Defect states and spin-orbital physics in doped vanadates: Yi.^Ca^VOs 
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We present a model for typical charged defects in weakly doped Yi-ajCa^VOa perovskites and 
study how they influence the magnetic and orbital order. Starting from a multiband Hubbard 
model we show that the charge carriers introduced by doping are bound to the Ca defects with large 
binding energy of ~ 1 eV at small doping, and give rise to the in-gap absorption band observed 
in the optical spectroscopy. The central position of a generic Ca defect with eight equidistant 
vanadium neighbors implies a partly filled defect band and permits activated transport due to 
Coulomb disorder. We explore the effect of bound charge carriers on the dynamics of the {yz,zx} 
orbital and spin degrees of freedom in the context of a microscopic strong-coupling model. After 
deriving the superexchange interactions around the doped hole we show that the transition from 
G-type to C-type antiferromagnetic order is triggered by the kinetic energy of doped holes via the 
double exchange mechanism. The defect states lead to local modification of orbital correlations 
within ferromagnetic chains along the c axis — some of them contain hole defects while the charge- 
orbital coupling suppresses locally {yz, zx} orbital fluctuations in the others. Thereby Ca defects 
provide a physical mechanism for spin-orbital dimerization along the ferromagnetic bonds, suggesting 
that in the C-AF phase of weakly doped Yi-^Ca^VOa dimerization increases with doping. 

PACS numbers: 75.10.Jm, 71.10.Fd, 71.55.-i, 75.25.Dk 



I. EXPERIMENTAL MOTIVATION 

Recent experimental and theoretical investigations of 
transition metal oxides have revealed the interplay be- 
tween spin, charge and orbital degrees of freedom, lead- 
ing to dramatic changes of magnetic and transport 
properties^ The best known examples of this joint effect 
of several degrees of freedom are high-temperature su- 
perconductivity and colossal magnetoresistance (CMR). 
When orbital degrees of freedom contribute, as in the 
CMR manganites, they may strongly influence magnetic 
properties and also may play a key role for the charge 
transport^— Particularly in the manganites the strong 
coupling between orbital states and local lattice distor- 
tions plays a prominent role, and orbital order is stabi- 
lized to a large extent by the Jahn- Teller interactions. 5 
Under these circumstances quantum orbital fluctuations 
are suppressed and classical orbital order determines the 
spin structure and dynamics £ This is most clearly dis- 
played by the large difference between the temperatures 
of the orbital and magnetic phase transitions in LaMnC>3, 
being T 00 = 780 K and T N = 140 KZ 

A challenge for the theory of spin-orbital systems and 
a qualitatively different situation is encountered in the 
-RVO3 perovskites, where i?=Lu,Yb,- • • ,La. In these 
perovskites, i.e., controlled by ti g valence electrons, the 
orbital-lattice coupling is weak while the spin-orbital cou- 
pling is the dominant interaction and thus orbital fluctu- 
ations are not suppressed. A common feature is the onset 
of the G-type alternating orbital (G-AO) order below the 
characteristic orbital ordering temperature Too which is 



in these compounds relatively low, Too — 200 K, and 
comes close to the Neel temperature 100 < Tni < 140 
K for the magnetic transition to the antiferromagnetic 
(AF) phase with AF order in ab planes accompanied by 
ferromagnetic (FM) order along the c axis, the G-AF 
phased Unlike e g orbitals in the CMR manganites, in 
this class of compounds the t2 g orbitals may form orbital 
ordered states which are subject to strong orbital quan- 
tum fluctuations. The consequences of quantum spin- 
orbital interplay in the -RVO3 perovskites were discussed 
in the theory^— and have been observed in several 
experiments.— ~— The changes of spin and orbital corre- 
lations are responsible for the temperature dependence of 
the optical spectral weights; 12 ' 15 as well as for the dimer- 
ization of FM interactions in the G-AF phase observed in 
the neutron scattering in YV03^ and also in LaV03i21 
Related instability of FM chains toward dimerization oc- 
curs at finite temperature when spin and orbital degrees 
of freedom couple. 22 Orbital fluctuations and their com- 
petion with orbital-lattice coupling play also a crucial role 
for the explanation of the nonmonotonous dependence of 
the orbital transition temperature on the radius of R 
ions along the -RVO3 series*^ 

The phase diagram of the perovskite-type -RVO3 
compounds 8 shows several spin- and/or orbital ordered 
phases. In the regime of compounds with low values of 
ionic radii t\r of rare earth ions R, another AF phase with 
complementary G-type AF (G-AF) order- 3 accompanied 
by G-type alternating orbital (G-AO) order (with stag- 
gered orbitals in ab planes and repeated orbitals along 
the c axis) appears below the second magnetic transi- 
tion at T N2 ^ for example in YV0 3 T N2 = 77 KMrZL In 
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addition, recent Raman experiments^ suggest that the 
short-range orbital fluctuations of the G-type occur in 
this intermediate C-AF phase — they coexist with the 
C-AO order and make it thus quite different from the 
one observed in LaVC>3. The magnetic exchange con- 
stants which determine the magnons in the C-AF phase 
are about a factor two smaller than those found in the 
low-temperature G-AF phased Therefore, it has been 
argued that this phase transition in YVO3 follows from 
large entropy of spin and orbital excitations in the G-AF 
phase^^ but the observed reduction of the energy scales 
of magnetic excitations remained puzzling and could not 
be explained by theory so farJ^ 

In recent years the effect of doping in several cu- 
bic vanadium oxides systems such as Lai-^Sr^VOa, 
Pr 1 _ x Ca a V0 3) Ndi-aSr^VOa, and Yi_ x Ca x V0 3 has 
been systematically explored by various experimental 
techniques ^tM^ Resistivity data, specific heat and mag- 
netic measurements have been used to set up the phase 
diagram as function of doping. 32 In contrast to the 
high-T c cuprates, where the metal-insulator (MI) tran- 
sition is found at a few percent doping ^ in the vana- 
dates the MI transition occurs at much higher doping 
concentrations 18% Sr in La^Sr^VOs and 
even up to 50 % Ca in Yi_ x Ca x V03. The evolution 
of optical spectra with doping for these two systems 
shows that the defects lead to impurity states which 
appear as absorption band deep inside the Mott gap. 33 
This suggests that bound small polarons are the cause of 
the MI transition at such high doping concentrations. 33 
It is eventually the growth of the mid-infrared absorp- 
tion with increasing doping and the gradual shift of this 
absorption toward zero energy which accompanies the 
insulator-metal transition. 

When taking all these experimental features into ac- 
count one arrives at a clear physical pictured (i) most 
importantly, the trends of the optical conductivity show 
that the edge of the Mott-Hubbard gap is essentially un- 
affected by doping, and only fades away when the MI 
transition is approached; (ii) defects play a central role, 
not just by introducing holes, but as generators of deep 
impurity states which appear in the optical conductiv- 
ity as midgap absorption at low doping; and finally (iii) 
the defects introduce two distinct energy scales. On one 
hand one finds the defect binding energy of about 1 eV, 
i.e., in the dilute doping regime, and on the other hand 
there is clearly an activation energy of ~ 0.1 eV or less in 
transport experiments i 34 ' 36 Our aim here is to show how 
generic Ca defects doped into the Mott-insulator YVO3 
explain these phenomena in a natural way. 

We begin our investigation with a discussion of the 
properties of Ca defects inserted into an orbital degen- 
erate Mott-Hubbard insulator. We adopt a multiband 
Hubbard model description of the ti g electrons and 
we use the unrestricted Hartree-Fock (HF) method i 39 i 
This approach is chosen as it allows us to introduce in a 
straightforward manner the lower Hubbard band (LHB) 
of the tig orbital states as well as the upper Hubbard 



bands (UHBs) with the appropriate multiplet splitting. 
The most pronounced effect of the defect is the Coulomb 
potential of the Ca-impurity which gives rise to an up- 
ward shift of the vanadium ti g states in the neighbor- 
hood of the defect. This leads to deep impurity states 
in the Mott-Hubbard gap. Interestingly, on one hand 
a Ca defect introduces one hole, but on the other hand 
it generates defect states on eight equivalent vanadium 
neighbors* 41 - Thus the topmost defect states that are split 
off the LHB are partially filled and pin the chemical po- 
tential. 

Transport inside the narrow defect band at weak dop- 
ing will be affected by the Coulomb disorder— of the 
charged defects. Consequently one expects Anderson lo- 
calization as was actually already conjectured by Mott 37 
for the Lai-^Sr^VOs system. Subsequently we proceed 
to our central aim, namely the investigation of the effect 
of defects on the spin-orbital dynamics in the dilute limit. 
Hence important aspects concerning defects, namely: (i) 
consequences of disorder, (ii) the role played by the long- 
range Coulomb interaction, and most importantly (iii) 
the interaction effects at higher doping, will not be dis- 
cussed here, but will be addressed elsewhere. 

One of the striking differences between the G-AF and 
G-AF phase is the stability of the latter phase in doped 
Lai-^Sr^VOa and Yi-^Ca^VOs compounds. For ex- 
ample, in Lai-^Sr^VOs the G-AF order survives even 
beyond the insulator-to-metal transition at doping x ~ 
0.18, 31 and disappears only at x ~ 0.2632 3 The G-AF 
of YVO3, however, is fragile and is destabilized in the 
Yi-^Ca^VOa compounds already at x ~ 0.02, where 
the G-AF phase takes overi 32 i 33 We shall argue below 
that in order to understand this behavior it is crucial to 
treat explicitly the dynamics of orbital degrees of free- 
dom. The cubic symmetry is broken at V 3+ ions due to 
the orthorhombic lattice distortion which occurs below 
the structural transition and favors energetically the elec- 
tron occupancy of xy orbita h 24 ' 27 This symmetry break- 
ing was also confirmed by the electronic structure cal- 
culations performed for LaV03 and YYO^.^— Due to 
Hund's exchange both electrons in a d 2 configuration at 
V ion form a high-spin S = 1 state, so the second ti g 
electron occupies either yz or xz orbital, resulting in a 
xy 1 (yz/zx) 1 local configuration at each V 3+ site. 

The orbital state is quite different in both magnetic 
phases of YVO3, below and above Tn2- Lattice distor- 
tions are large in the low-temperature G-AF phase and 
suggest G-AO order. This orbital order is further sta- 
bilized by increasing pressure^ 7 - Above Tn2 the distor- 
tions decrease and are compatible with a weak G-type 
AO (G-AO) order £1 It was suggested by Ishihara 4 ^ that 
the phase transition at Tn2 could originate from the or- 
bital physics and would be triggered by orbiton softening 
induced by the reduction of the spin order parameter. 
While the orbital degrees of freedom certainly play a role 
as the orbital order indeed changes at T/V2, there is no 
evidence of orbiton softening so far. It could be expected 
that the observed transition is caused instead by local 
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FIG. 1: (Color online) Schematic view of the lattice of vana- 
dium sites with occupied tig orbitals and a single hole h in- 
troduced by doping a Ca defect D (sphere in the center) in 
the C-type orbital structure of the low-temperature phase in 
weakly doped Yi-^Ca^VOs. For clarity the Y and O ions and 
the orbital phases are not displayed. The hole occupies pref- 
erentially one of the V sites that form a cube around the Ca 
impurity due to the attractive Coulomb potential of the im- 
purity. At the undoped V sites only the topmost occupied t2 g 
orbitals are shown. The c) = \xy) orbitals occupied at each 
V ion and residing at lower energy are not shown. Orbital 
polarization, Eq. (|3.12|l . distorts the \a) = \yz) (|6) = \xz)) 
C-type alternating orbital order in the neighborhood of the 
defect (dashed box) and favors occupation of one of the two 
{|+)»j |— )i} rotated orbital states that minimize the orbital- 
defect interaction (see Sec. IIII D|) . 



phenomena close to Ca impurities in Yi_ x Ca x V03 rather 
than by the global change of orbital excitation scale. We 
suggest that it is plausible that impurities could locally 
destabilize the C-AO order, and introduce a microscopic 
model to treat this effect below. An earlier theoretical 
analysis within the charge-transfer model has shown that 
both phases are indeed energetically close,— and one may 
thus expect that small changes of the thermodynamic po- 
tential around Tjv2 could trigger a first order magnetic 
phase transition. 

The purpose of this paper is to investigate the local 
changes in the electronic structure caused by Ca defects 
in the sublattice of Y ions, and to derive the micro- 
scopic model leading to a scenario which could explain 
why the G-AF phase is so fragile under Ca doping in 
Yi_ x Ca x V03. Optical spectroscopy has shown that the 
absorption in the Mott-Hubbard gap grows as a function 
of Ca doping in Yi_ x Ca x V03^^ In the low doping 
range x < 0.05 it has a characteristic energy close to 
1.2 eV,— and moves to lower energies at higher doping. 
The intensity of the mid-infrared peak increases propor- 



tionally to doping x. This new excitation energy which 
occurs in the doped system provides a measure of the 
binding energy of a doped hole which confines the hole 
to the V§ cube in the immediate neighborhood of the 
charge defect, see Fig. [TJ For the undoped bonds of this 
cube we invoke the mechanism of local polarization of 
orbitals near the Ca 2+ charge defects in the Y 3+ back- 
ground. It will be shown that such defects, together with 
holes introduced into the ti g orbitals of vanadium ions, 
trigger the collapse of the G-AF order. 

The major part of our study is based on an effective 
low-energy Hamiltonian which can be denoted as a spin- 
orbital t-J model; the model provides a description of 
magnetism and hole motion in the lower Hubbard band of 
a doped and orbital degenerate Mott-Hubbard insulator. 
In the undoped case the model coincides with the spin- 
orbital superexchange model^ ' 10 ' 12 " — which was shown 
earlier to lead to a satisfactory description of spin-orbital 
physics, e.g., in YVO3. The spin-orbital t-J model pro- 
vides moreover the description of the kinetic energy of 
doped holes and describes the effect of charge impuri- 
ties on the holes and the spin-orbital degrees of freedom. 
In the concrete case of Yi-^Ca^VOa we show that the 
changes introduced by the defects have severe effects on 
the orbital dynamics, namely: (i) are responsible for the 
destabilization of the coexisting G-AF and G-AO order 
with increasing doping x, and (ii) lead to an enhanced 
tendency towards dimerization in the weakly doped G- 
AF phase. 

The paper is organized as follows. First, in Sec. m 
we describe Ca 2+ charge defects in Yi-^Ca^VOa within 
the frame of a degenerate Hubbard model for t^ g orbitals 
and analyze the effect of defects on the densities of states 
derived within the HF approximation. In the follow- 
ing Sec. IIIII the spin-orbital t-J model is outlined. The 
model Hamiltonian contains apart from spin-orbital su- 
perexchange interactions also pure orbital interactions in- 
duced by the lattice, see Sec. IIII Bl The effective double- 
exchange mechanism is introduced in Sec. IIII CI Next we 
derive the modified superexchange in the vicinity of the 
hole in Sec. IIII El and develop the idea of orbital polariza- 
tion around charge defects in Yi^-jCa^VC^ in Sec. IIII Dl 
With these terms, which supplement the spin-orbital 
model for undoped cubic vanadates, the complete t-J 
model poses a rather complex many-body problem, and 
we derive from it two effective one-dimensional (ID) em- 
bedded orbital models for the G-AF and G-AF phase sep- 
arately, see Sees. IIVBI and IIVCI Thereby we treat spin 
correlations of S — 1 spins in a classical approximation 
and focus on the orbital dynamics that plays a decisive 
role for the observed phase transition in Yi-^Ca^VOs. 
Here we show how the orbital polarization (see also Ap- 
pendix A) develops around charge defects in the orbital 
chains for both magnetic phases. Next we consider a hole 
in both magnetic phases (Sec. [Vj) , and derive the ID or- 
bital chain models appropriate for the G-AF phase and 
G-AF phase in Sees. IV Al and IV Bl Interactions around 
the hole are described in the G-AF phase by the d 2 — d 1 
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superexchange (Appendix B) and the kinetic energy is 
controlled as in a typical double exchange mechanism 
by the spin orientations (discussed also in Appendix C). 
This implies, e.g., free hole hopping along the c axis in the 
C-AF phase due to the FM alignment of spins. In Sec. 
IVII we introduce a statistical treatment of orbital corre- 
lations and energy contributions at finite doping (Sec. 
IVI A[) and come to the conclusion that defects trigger 
dimerization of orbital correlations and of FM spin ex- 
change constants along the c axis in the C-AF phase 
(Sec. IVIBI) . Finally, using the effective models derived 
for both magnetic phases we present a scenario which ex- 
plains why the phase transition from the G-AF to dimer- 
ized C-AF phase takes place already at small doping, see 
Sec. IVI CI A summary and outlook are given in Sec. I VIII 

II. DEFECT STATES IN RVO :i PEROVSKITES 

A. Degenerate Hubbard model for t2 g electrons 

We begin with analyzing the consequences of strong 
Coulomb interactions in the framework of a multi-orbital 
Hubbard model introduced here to describe the doped 
Yi-^Ca^VOs compounds. This model will later form 
the basis for the derivation of an effective superexchange 
and a related spin-orbital t-J model for the orbital de- 
generate case, see Sec. MI Al In the ionic picture of the 
undoped YVO3, V 3+ ions are in a d 2 electronic configu- 
ration with partly filled t2 g orbitals. As we show below, 
the qualitative features of the optical spectra, i.e., the 
multiplet splitting of the Hubbard bands and the posi- 
tion of defect states, may be reproduced by considering 
a multiband Hubbard model for ti g electrons, 

H = H t + Hqf + Hi n t + -ffimp , (2.1) 

where the four terms stand for the kinetic energy {H t ), 
crystal- field (CF) splitting (Hqf), local Coulomb inter- 
actions (flint), and Coulomb potential generated by a 
Ca impurity (fli mp ). On one hand this model can also 
serve as a basis for the qualitative discussion of the pho- 
toemission (PES) and inverse PES in the vicinity of the 
Mott-Hubbard gap. On the other hand, it provides the 
basis for the derivation of the spin-orbital Hamiltonian 
which serves for a transparent description of the magnetic 
and orbital structure, as well as of spin and orbital exci- 
tations. Our aim here is to explore further the changes 
of the excitation spectra resulting from the presence of 
Ca 2+ defects in the lattice of Y 3+ ions and their impact 
on the vanadium ti g electrons, with help of this simplified 
Hamiltonian. 

The kinetic energy is given by: 

7 <ij>||7.«(7).CT 

(2.2) 

where d\ aa is electron creation operator for an electron at 
site i in orbital state a with spin a =t, i- The summation 



runs over three cubic axes, 7 = a, b, c, the bonds (ij)ll7, 
and the hopping t a conserves the t2 g orbital flavor. The 
effective hopping t a originates from two subsequent d — p 
hopping processes via the intermediate 2p v oxygen or- 
bital along each V-O-V bond. In principle it can be de- 
rived from the charge-transfer model with p—d hybridiza- 
tion t p d and charge-transfer energy and one expects 
in the present case t — t 2 d / A ~ 0.2 eVM- Only two ti g or- 
bitals, labelled by a(j), are active along each bond 
and contribute to the kinetic energy Eq. (|2.2[) . while the 
third one lies in the plane perpendicular to the 7 axis 
and the hopping via the intermediate oxygen 2p v oxygen 
is forbidden by symmetry. This motivates a convenient 
notation used below, 

\a) = \yz), \b) = \xz), \c) = \xy), (2.3) 

where the orbital inactive along a cubic direction 7 is 
labelled by its index as I7). 

In agreement with a commonly accepted picture^ 
information obtained from the electronic structure 
calculations ) 44 " 46 and with the results obtained using the 
point charge model^ we assume that the xy (c) orbitals 
are energetically favored and thus occupied and inactive 
at low temperature, while the remaining yz and xz or- 
bitals are nearly degenerate and represent the t2 g orbital 
doublet, with both orbitals active for the hopping (and 
the superexchange) along the c cubic axis. 

The nonequivalence of the t2 g orbital states is de- 
scribed by a CF splitting term which favors the c orbitals, 

flcF = -^ A «"«^' ( 2 - 4 ) 

where n, Q(T = d\ aa d iaa is an electron density operator, 
and A Q = (5 QiC A c with A c > . As a result, when 
A c > t as we estimated^ c orbitals are filled by one 
electron at each site in a strongly correlated system, and 
the second electron occupies one of the orbitals in the 
{a, b} doublet, leading to the c\(a,b)\ configuration at 
each site i. This broken symmetry situation corresponds 
to electron densities 

n ic ~ 1, n ia + n ib ~ 1, (2.5) 

within t2 g orbitals at undoped V 3+ ions. Note that one 
expects that the cubic symmetry with n^ 7 = 2/3 is re- 
stored at high temperature, but this situation will not be 
analyzed here as it has no influence on the mechanism 
of the phase transition from the G-AF to C-AF phase 
which occurs in Yi-^Ca^VOs at low temperature under 
increasing doping. 

The local Coulomb interactions between t2 g electrons 
at V 3+ ions are described by the degenerate Hubbard 
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Hamiltonian^i with the interacting part: 

flint = U ^ Uia-fTliai + ( U — -Jh J ^ «ia™i/3 
ia ^ ' i,a</3 

i,a</3 

- 2J# 2J ^ia ' Si/3 . (2.6) 
i,a</3 

Here rii a — rii aa is the corresponding electron density 
operator in orbital a at site i, and spin operators Si a — 
{Sf a , Sf a , Sf a } are related to fermion operators in the 
standard way, i.e., 

^ia — diafdial > Sf a = — (riia^ — rijaj.) . (2-7) 

The first term in Eq. (|2.6j) describes the intraorbital 
Coulomb interaction U between electrons with antipar- 
allel spins. The second term stands for the interorbital 
Coulomb (density) interaction, the third one is called 
frequently the "pair- hopping" term, and the last one is 
Hund's exchange Jh- The choice of coefficients in Eq. 
(|2.6p guarantees that the interactions satisfy the rota- 
tional invariance in the orbital space. 51 This Hamilto- 
nian is exact when it describes only one type of 3d or- 
bitals which are partly occupied, as ti g orbitals in the 
present case of the -RVO3 perovskites, and the inter- 
actions are then given by two parameters: (?) the in- 
traorbital Coulomb element U and (ii) the interorbital 
(Hund's) exchange element Jh- These elements may be 
expressed by the Racah parameters {A, B,C}. For ti g 
electrons one finds £>£2 

U = A + AB + ZC, (2.8) 
J H = 3B + C. (2.9) 

Finally we introduce the Coulomb interaction between 
the t2 g electrons at a V-site r, and the effectively negative 
charged defects — with charge Qd = e at site R n: 

iec„ 

where rii = J2 a n ia is the total ti g electron density. Here 
i (z Ci (i £ Coo) denotes a calculation where in the sum 
only nearest (all) V neighbors of a defect are considered. 
The Coulomb potential itself is long-ranged, 

W(r) = —, (2.11) 

and screened by the dielectric function e c of core elec- 
trons. We identify e c with the high frequency dielectric 
constant which, e.g., for YVO3 lies in the interval e(w) — 
5.0 ± 0.3 in the whole frequency range 0.8 < ui < 3.5 
e y54 rjv^g mog £ pronounced effect of this potential term 
is an upward shift of the vanadium states in the neigh- 
borhood of the defects. This gives rise to bound states 



in the Mott-Hubbard gap. It is important to realize that 
on one hand a Ca-defect D introduces one hole, but on 
the other hand it generates defect states on eight equiva- 
lent vanadium neighbors. Thus the topmost defect states 
that are split from the LHB are only partially filled and 
pin the chemical potential [i. 

In general the coordinates of the defects R n will be sta- 
tistically distributed. In fact, defects will also feel some 
repulsion and avoid clustering. We will not explore these 
aspects here, as we are concerned mainly with the dilute 
doping regime. Nevertheless it is clear that transport is 
strongly affected by disorder and results from the motion 
of holes in the defect band. 

Furthermore, in the immediate neighborhood of a Ca 
defect the strength of Coulomb interaction influences the 
orientation of t^g orbitals filled by electrons. This effect 
modifies the orbital state and will be described below via 
an extra crystal field term which acts on the orbitals of 
the V ions in the vicinity of a defect (see Sec. IIII D[) . Such 
more subtle effects concern the low-energy behavior and 
are analyzed in the framework of the effective spin-orbital 
t-J model in Sec. IIIII Therefore, we do not analyze it 
here, similar as the intersite orbital interactions which 
originate from the distortions of VOe octahedra. 



B. Unrestricted Hartree-Fock analysis of the 
degenerate Hubbard model 

Qualitative insight into the correlated electronic struc- 
ture of partly filled t 2g states can be obtained from the 
HF approximation. When the HF approximation is used, 
the " pair- hopping" term does not contribute and the lo- 
cal spin exchange interactions S{ a ■ Sip contribute only 
with the Ising term cx SfSjf& i.e., one may use an ap- 
proximate expression, 

— 2S ia ■ Sip ~ - V] (nia^n^g- — nio^riipa) . (2.12) 



As a result, local electron interactions Eq. (|2.6I) are given 
approximately by electron density operators {rii aa }: 

H int ~ 

+ (U-2J H ) Y, 

ia i,a<{3,a 

+ (U-3Jh) Y, n ia<7mp<T- (2-13) 

i,a<@,a 

The form given in Eq. (|2.13[) violates the rotational in- 
variance of local Coulomb interactions, 51 but is sufficient 
for our purpose since all the terms which appear in the 
HF approximation arise from iti^ 5 - One finds the following 
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effective one-particle problem in an effective field, 

+ (U-2J H ) ^ [ni a a{nips) + {n ia a)nip s ) 

i,cx<Cf3,<T 

+ (U - 3J H ) £ 

i,ce<C{3,<r 

(2.14) 

The HF potentials are determined by the average densi- 
ties {(n iacr )} when H mt is replaced by H^f in Eq. (|2.13l) . 

We emphasize that although the quantum effects such 
as spin fluctuations and "pair hopping" of double oc- 
cupancies are neglected, the essential features of the 
Coulomb interaction are reproduced. This can be seen 
by considering charge excitations — > c^ n+1 ' ) c^ n-1 ' ) 
along a given bond (ij) in the lattice of transition metal 
ions with d n electronic configuration. To analyze excited 
states we assume that the electron number n is lower 
than the half- filled shell, i.e., n < 5 for the 3d shell (be- 
low we focus on n = 2 for the present problem of V 3+ 
ions in YVO3, where the ti g shell is half- filled at n = 3). 
There are three types of S n+1 ^ excited states: (i) high- 
spin (HS) states with all electrons in the same spin state 
(realized for a FM (ij) bond); (ii) low-spin (LS) states 
with all orbitals being singly occupied, and (iii) LS states 
with one doubly occupied orbital. The energies of these 
excitations are: 

£hs = U-3J H , (2.15) 
E$ = U - 2 J H + (n - 1)J H , (2.16) 
Ej$ = U+(n-l)J H ., (2.17) 

The HS excitation energy obtained in the HF approxi- 
mation reproduces the exact value obtained by the diag- 
onalization of the atomic ion Hamiltonian Eq. (|2 . 6[) . The 
remaining energies are systematically lower by Jh than 
the exact values 

4s,ex = U-J H + (n-l)J H , (2.18) 
^Sex = U + J H + (n-l)J H , (2.19) 

as the quantum fluctuation effects (spin- flips and "pair- 
hoping") were neglected. 

In the relevant regime of parameters for the -RVO3 per- 
ovskites, the electrons in the undoped YVO3 are local- 
ized in a Mott insulator. We consider here representative 
parameters^ with [7 = 4 eV and Jh = 0.6 eV. In this 
parameter range one may simplify the problem of calcu- 
lating the electronic structure in the HF approximation 
as the c orbitals are occupied at each site by one t2 g 
electron, and the magnetic state is determined by their 
magnetic moments. This follows from large Hund's ex- 
change which dictates that the spin direction of both t2 g 
electrons (c electron and (a, b) one) agree at each site. As 
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FIG. 2: (Color online) Partial density of states N(u>) (in eV - 
units) for t2 g states with {yz, zx} symmetry of the Mott in- 
sulator YVO3 obtained by the unrestricted HF calculation. 
The multiplets represent the occupied lower Hubbard band 
(LHB) while the unoccupied upper Hubbard band consists 
of high-spin (HS) and low-spin (LS) multiplet states, respec- 
tively. The Mott-Hubbard gap Amh separates the LHB and 
the HS band. The dashed step-like curve indicates the inte- 
grated number of states f(io) Eq. (|2.22p normalized to one. 
Parameters: U — 4, Jh ~ 0.6, t c — t a — tt — 0.2 (all in eV). 



a result, one obtains the high-spin S = 1 state at each 
V 3+ ion. 

We discuss first the numerical calculations of the den- 
sity of states N(u) for the undoped system. The results 
were obtained by solving the equations derived using the 
HF approximation, see Eq. (I2.14j) . on a cluster with 
Nq = L x x L y x L z sites (here we use L x = L y = L z = 10) . 
We considered the c orbitals singly occupied and ordered 
as in the C-AF phase of the undoped YVO3. The second 
electron at each site occupies one of the remaining ti g or- 
bitals {a, b}. The partial density of states for the {a, b} 
orbital doublet with eigenenergies e m in the cluster, 

%)^^-e ra ), (2.20) 

m 

is shown in Fig. [5] The chemical potential /1 in the 
undoped case is naturally chosen in the middle of the gap: 
fi = (ejv + e JV +i)/2, where ejv is the highest occupied 
energy for Nq electrons in {a, 6} orbitals in the system. 

The spectra are characterized by four subbands: (i) 
the LHB centered at £"lhb and (ii) the UHB which itself 
is split into three subbands (multiplet structure) corre- 
sponding to the HS excitations at -Ehs = U — 3Jh Eq. 
(|2.15[) , and two LS transitions centered at E^s = U — Jh 
Eq. (|2~T6|) and E$ = U + J H Eq. (|2T6]) . respectively. 
These energies are relative to £"lhb ■ HF results for these 
excitations calculated for U = 4, Jh = 0.6 and to = 0.2 
eV (t a = tb = t c = to) are listed in Table I and compared 
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TABLE I: Excitation energies of YVO3 inferred from optical 
dat a 33 ' 57 and calculated from the present theory. Here Amh 



is the Mott-Hubbard gap, and Ehs, and E^g are the 

energies (all in eV) of the high-spin and the two low-spin 
transitions measured from the center of the lower Hubbard 
band £?lhb- The last column gives the defect energy En 
relative to -Elhb at x = 0.02 Ca-doping. 
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with experimental values for YVO3 obtained by Fujioka 
et al~ and by Kupersbusch 5 £ deduced from optical spec- 
troscopy and ellipsometry, respectively. Interestingly the 
LS-transitions found in the ellipsometry study find a rea- 
sonable correspondence within the HF calculation. These 
values are also consistent with the underlying multiplet 
splitting of YVO3 in the optical spectral weight study of 
Ref. EI 

Next we turn to the Mott-Hubbard gap (see Fig. ^ 
which may be expressed as: 

Aun^U -3J H -W cSl (2.21) 

where W c s is the effective bandwidth of a Hubbard 
subband. A significant reduction of W e ff compared to 
the free bandwidth W is well known from the single 
band Hubbard model^ In a recent study a LHB width 
Weff « 3to was estimated for the incoherent hole motion 
of tig electrons^ Thus with the choice to — 0.2 eV we 
obtain for the Mott-Hubbard gap Amh — 1.6. We also 
note, that our estimate of the bandwidth of the LHB in 
the HF calculation is much smaller, namely Wjj ~ to- 
This is due to the neglect in the HF approximation of 
processes describing the incoherent motion of holes. A 
related interesting quantity that can be inferred from 
optical spectroscopy is the width of the HS-transition, 33 
r HS - 1.3 ± 0.2 eV in YV0 3 . As the width of the optical 
transition should be determined by the convolution of the 
LHB and HS-band, one expects Ths — 2W c ff. Thus Amh 
and the width Tjjs nn d a natural explanation in terms of 
the multiplet splitting and the effective Hubbard band- 
width Weff. 

The counting of states that contribute to the partial 
density of states N(lu) is straightforward; there is one 
electron per site which can occupy either a or b orbital 
and there are two spin flavors. Thus the filling is /(/tt) = 
1/4 as the total number of states is N t — 4iVo, where No 
is the number of sites, and there are No electrons that 
occupy the LHB. All subbands of the UHB have also 
the same weight of No each. This can be seen from the 
integrated and normalized density of states 

+UJ +00 

/(w)= J du'N{uj')/ J du)'N(u'), (2.22) 

—00 —00 



FIG. 3: (Color online) Spectral weight of Hubbard subbands 
per site at hole doping x for: (a) the single band Hubbard 
model and (c) Hubbard bands of doubly degenerate model 
with two spin and two orbital flavors. At doping concen- 
tration x the number of states (1 — x) in each subband of 
the UHB is determined by the number of filled states in the 
LHB. Whereas in (a) a hole can be filled either by an up- or 
a down-spin electron, as shown for a representative electron 
configuration in (b), in the orbital degenerate case (c) there 
are 4 choices to fill up the hole in the atomic limit, i.e., 4x 
empty states in the LHB for a system with x doped holes. 



which is displayed as the dashed curve in Fig. [2] While 
we use a broadening parameter 7 = 0.1 eV to smoothen 
the density of states N(w), no such broadening is used in 
the calculation of f(u)). 

It is well known that the number of states in the UHB 
depends on the hole doping ) 58 i 60 i.e., varies with electron 
filling in striking contrast to bands of ordinary semicon- 
ductors or insulators. It has also been realized that by 
creating a hole in a Mott insulator (with no orbital de- 
generacy) actually two empty states are created in the 
LHB, see Fig. (3Ja) as there are two options to fill up the 
holei^ir— That is, at a concentration of x holes, there 
is a fraction of 2x unoccupied states in the LHB above 
the chemical potential, and simultaneously the number 
of states in the UHB is reduced by a factor 1 — x. Such 
a spectral weight transfer upon electron doping has been 
recently observed in TiOCli^ 

In case of the present model with two orbital flavors 
{a, b} in addition to spin, there is one LHB and three 
UHBs, all of them with the same weight one at x = 0. 
At finite doping x > each added hole moves an empty 
state above the Fermi energy in the LHB and generates 
also three other unoccupied states in the LHB which are 
taken from the subbands of the UHB, see FigJ3fc). This 
corresponds to four possibilities to fill up a hole (with spin 
and orbital flavor) and thus 4x empty states belong to the 
LHB at doping x. Yet it is also important to emphasize 
here that this does not imply that there are really Ax 
free states in the LHB that can all be simultaneously 
occupied. Actually there is only space for x electrons, as 
with each electron added three states are shifted back to 
the UHBs. 
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FIG. 4: (Color online) Partial density of states of a doped 
Mott insulator with degenerate ti g orbitals {a, 6} at x — 0.02 
doping concentration. Partially filled defect states D are split 
off from the LHB by the defect potential Vd = 1 eV. The 
vertical dashed line indicates the filling of the system at fi- 
nite doping; D* marks defect states split from the highest LS 
excitation; further notations and parameters as in Fig. 2. 



C. Hartree-Fock calculation for defect states 

The major change in the spectra induced by doping is 
the occurrence of the defect states D inside the Mott- 
Hubbard gap. Figure [4] shows the density of states at 
2% Ca-doping. The spectra were obtained by a calcula- 
tion using the HF approximation for well separated de- 
fects acting on ti g electrons with the defect potential, Eq. 
(|2.10j) , which includes the interaction with the eight near- 
est neighbor V ions at distance d\ = \/Z dyo (dvo — 2.0 
A is the vanadium-oxygen distance). The estimate of 
the nearest neighbor defect potential Vd appropriate for 
YV0 3 is straightforward, i.e., W{d x ) = Vb ~ 1.0 eV, 
where we used the dielectric constant of core electrons 
e c ~ 5.0. The summation in Eq. (|2.10l) was restricted 
to i G Ci, thus only a short-range defect potential was 
included, and disorder effects do not come into play. 
Thereby the interpretation of spectra and energy scales 
is simplified and more transparent for our purposes. 

The calculations show that for large Vd ~ 1.0 eV the 
energy Ed of the defect states D relative to the center of 
the LHB is approximately given by Ed ~ Vd (see Table 
I). The complementary splitting of the defect states from 
the center of the HS-band is E' D ~ 1.2 eV, and appears 
also consistent with optical spectroscopy data<22 Depend- 
ing on the value of the defect potential, the defect states 
D appear either at the upper edge of the LHB, or develop 
to a separated maximum within the Mott-Hubbard gap, 
see Fig. [SJ Figure GDJb) nicely shows that each UHB 
(HS, LS (1) and LS (2) ) has its own defect satellite. Such 
states can be observed, however, only when a satellite of 




FIG. 5: (Color online) Partial density of states N(u}) of a ti g 
Mott insulator at x = 0.02 doping concentration for different 
values for the defect potential: (a) Vd = 0.1 eV, (b) Vd = 0.5 
eV, and (c) Vd = 1.0 eV. The vertical dashed line indicates 
the chemical potential n which is pinned by the defect states 
D. The binding energy Eb of the defect states measured from 
the edge of the LHB grows with increasing Vd ■ An alternative 
measure is the defect energy Ed (E' d ) relative to the center 
of the LHB (HS-band), respectively. Parameters: U = 4, 
Jh = 0.6, t c =t a =t b = 0.2 (all in eV). 



a given Hubbard subband is well separated from the next 
subband. 

We emphasize that in the model with two orbital fla- 
vors each Ca 2+ defect introduces one hole into eight de- 
fect states that are split off from the LHB. The filling 
fraction at doping concentration x — N x /Nq is therefore, 

, , N - N x 1 , 
f(u) = — - = -(1-0:), 2.23 

which fixes the chemical potential /x. The total number 
of occupied electron states is now simply obtained by 
multiplication with 47Vq: 



N occ = (1 - x)N , 



(2.24) 



and the number of holes is N h — xN . The number 
of states in the LHB, however, is reduced to JVlhb = 
(1 — 8x)No due to the appearance of Nd — 8xNq defect 
states, of which Nff c — 7xNo are occupied. 

Although there is a similarity to defect bands in doped 
semiconductors a striking difference is that the defect 
states here are completely derived from the LHB. Hence 
the defect band in Fig. [5jc) would be fully occupied if no 
holes were added. However each added defect generates 
not only the defect states but also contributes a single 
hole. We suggest that it is the motion of the holes in 
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the defect band, i.e., via hopping from one defect to the 
next one which occurs in presence of Coulomb disorder 
that explains the small excitation energies observed in 
transport experiments ! 34 ' 36 



ates an attractive potential and confines a hole to the 
immediate neighborhood of the charge defect. 

Strong electron correlations in the Yi-^Ca^VOs com- 
pounds motivate the spin-orbital t-J model, 



III. SPIN-ORBITAL t-J MODEL 
A. Model Hamiltonian for Yi-^Ca^VOa 

We now turn to the derivation of an effective low- 
energy Hamiltonian which describes the interactions of 
the spin and orbital degrees of freedom as well as the 
motion of doped holes. For the undoped case the effec- 
tive microscopic spin-orbital model that describes elec- 
trons in t2 g orbitals of V ions has already been worked 
out£ ' 10 i 12 ~ — In the undoped compound one deals with 
the t\ g electronic high-spin (5 = 1) state at each V 3+ 
ion, and an {a, b} orbital degree of freedom. 

When a correlated insulator with active orbital degrees 
of freedom is doped, rather complex processes occur when 
holes are doped. The motion of a hole may generate spin, 
orbital, or spin-orbital excitations on its path i 50 i 59 i 65 ~ — 
Therefore, designing a microscopic model which captures 
the essential physical mechanisms in such a situation and 
is mathematically tractable at the same time is a nontriv- 
ial and demanding task. Below we introduce such a mi- 
croscopic model to describe the changes of magnetic and 
orbital order that occur by hole doping in Yi_ x Ca x V(_3. 
It stems from the spin-orbital model for the undoped 
-RVO3 perovskites and includes the superexchange and 
the orbital interactions induced by the lattice^ii 3 . As in 
Sec. II, the kinetic energy arises from the hopping t be- 
tween two t2 g orbitals at neighboring V sites, and the or- 
bital flavor is conserved^ Electron-electron interactions 
are described by the degenerate Hubbard models Eq. 
(|2.6|) , written in this case for two t 2g orbitals {a, b}, with 
intraorbital Coulomb element U and Hund's exchange 
J H , see Eqs. flUBJ} and ([2~3]) . 

As in other Mott (or charge-transfer) insulators charac- 
terized by the realistic regime of parameters with t <C U, 6 
the magnetic and orbital interactions of strongly corre- 
lated t 2g electrons in the -RVO3 perovskites are then de- 
scribed within the spin-orbital superexchange Wj, with 
the superexchange constant 



A realistic model for the undoped -RVO3 perovskites con- 
tains also the orbital interactions "Horb which follow from 
the orbital-lattice coupling and are responsible both for 
the energetic proximity of the G-AF and C-AF phases 
in YVO3,— 1 10 ] 13 and for the systematic trends observed 
for the orbital and magnetic phase transition of the se- 
ries of -RVO3 perovskitesfi 4 - Here we present an appropri- 
ate extension of this model adequate for weakly doped 
Yi_j.Caj.VO3 compounds, which includes the charge- 
orbital interactions around the Ca 2+ defects. It gener- 



H t2g = V (Ht+HHund+Hj+Horh+Hi mp +HD+Hi) V , 

(3.2) 

where V are the projection operators which remove triply 
(and higher) occupied V ions and guarantee that either 
the {a, b} orbital doublet is singly occupied (at V 3+ ion) 
and spin is 5 = 1, or the {a, 6} doublet is empty (at 
V 4+ ion), i.e., it contains a hole generated by doping. 
The spin-orbital superexchange for d 2 — d 2 pairs of V 4+ - 
V 3+ ions (Hj), and the orbital-lattice interactions (H ib) 
stand for the effective strong coupling model^ see Sec. 
IIIIB1 that was used before to explain the temperature 
variation of optical spectra^- and the phase diagram 1 ^ 
of the i?VC>3 perovskites. When holes are doped, several 
other terms are needed: (i) the hopping of {a, b} elec- 
trons in the restricted space (% t ) Eq. (|2.2[) ; (ii) Hund's 
exchange "HHund between s = 1/2 spins of an {a, b} elec- 
tron and a c electron at site i; (iii) the impurity potential 
Himp given by Eq. (|2.10[) : (iv) the polarization interac- 
tion at the V sites near the charge defects {T~Ld)\ (v) the 
superexchange for d 2 — d 1 pairs of V 4+ -V 3+ ions (?_/). 
These terms are introduced in Sees. IIII CiljHDl and lHI El 
respectively. 



B. Superexchange in the undoped YVO3 

The third term in Eq. (|3.2p is the spin-orbital su- 
perexchange 'H.j. The anisotropic electron distribution 
between the ab planes and the c axis, see Eqs. (12.51) . 
is responsible for a particular form of the superexchange 
Hj, with broken cubic symmetry ^ 10 ' 12 In fact, this sym- 
metry breaking is responsible for strong {a, b} orbital 
fluctuations which stabilize the C-AF phase with FM in- 
teractions along the c axis. The superexchange interac- 
tions between two V 3+ ions in the undoped YVO3 with 
5 = 1 spins at sites i and j arise from virtual excitations 
dfd 2 =i d^dj along the concerned bond (ij), promoted by 
the hopping t which couples pairs of identical active t 2g 
orbitals. A single hopping process generates a d 3 con- 
figuration at site i, either with three different orbitals 
occupied by a single electron each, or with a double oc- 
cupancy in one of the two active orbitals. Therefore, the 
d\ excited state may be either a high-spin 4 A 2 state with 
energy (U — 3Jh), see Eq. (12.151) . or one of three low-spin 
states: 2 E, 2 T 1 or 2 T 2 with energies^ 6 - U and (U + 2J H ) 
— they are shown in Fig. 1 of Ref. H, see also Eqs. (|2.18[) 
and (pTT9l) . 

This perturbative consideration leads to the spin- 
orbital superexchange model for 5 = 1 spins in cubic 
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where (ij) is a nearest neighbor bond and J is the su- 
perexchange constant given in Eq. (13.11) . The superex- 
change Eq. (|3.3p follows from the degenerate Hubbard 
model Eq. (|2.ip as described in Ref. [l3| for the case when 
c orbitals are occupied, see Eq. (|2.5p . and a second elec- 
tron occupies the {a, 6} doublet at each site, i.e., n.; = 1 
with 



the G-AF order in the undoped YVO3. These interac- 
tions increase along the -RVO3 perovskites towards the 
compounds with small ionic radii, and it happens to be 
just for YVO3 that they tip the balance between the two 
types of magnetic order in favor of the G-AF phase which 
is more stable at low temperature^ 



C. Effective double exchange model 



The first term H t = VH t V in Eq. ([52} is the kinetic 
energy which after projection describes only hopping pro- 
cesses within the LHB. In the LHB only an electron at 
nearest neighbor site of the hole can hop by interchang- 
ing its position with the hole^i The remaining hopping 
processes describe either excitations to the UHBs that 
are included in the superexchange Hj between two V 3+ 
ions, or low spin charge excitations at the hole site that 
contribute to the superexchange Hi for V 4+ -V 3+ pairs 
of ions, see below. 



Til — ll'ia 1 Tlib- 

It depends on Hund's exchange, 

Jh 



(3.4) 



(3.5) 



due to the charge excitations to the states of the UHBs 
described above and shown in Fig. [31 via the coefficients 
(the coefficient 7*2 = 1 arises for the intermediate energy 
LS excitations at V +2 ions), 
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l + 2ry 



(3.6) 



The operators = {rf ,rf , rf} describe orbital r = 1/2 
pseudospins defined here (for each direction 7 = a, b, c) 
by the doublet of active {a, b} = {yz, xz} orbitals, and 
are given by the Pauli matrices, i.e., 



1 y 
2°*' 



(3.7) 



As both orbitals are active along the c axis, the orbital 
part is then given by a scalar product T{ ■ fj. 

The orbital-orbital interactions which follow from lat- 
tice distortions of both the Jahn- Teller type and GdFe03- 
type are of the {orm&±&mi 



Ho 



v ab 



(ij)\\ab 



Vc 



(y)ll c 



(3.8) 



The interactions in the ab planes V a b > follow from the 
Jahn- Teller distortions, while the ones along the c axis 
V c > favor the G-AO phase and thus help to stabilize 



Further restriction on the hopping is introduced by the 
breaking of cubic symmetry in doped Yi-xCa^VOs by 
the actual anisotropic electron distribution over the ti g 
orbitals given by Eqs. (|2.5p . This leads to the symmetry 
breaking between the bonds in ab planes and along the 
c axis. 12,13 As we have discussed in Sec. [Hi hole doping 
occurs in the orbital doublet {a, b}, and the c orbitals 
are filled also in doped systems by one electron each. 
Therefore, the c electrons are immobile in the strongly 
correlated regime and contribute only to virtual d^d 2 ^ 
d^dj excitations which generate the superexchange along 
the considered bond (ij) \\ ab. 

In the large U regime (U 3> t) the kinetic energy of 
the {a, b} electrons is finite only in a doped system when 
the hopping process may occur in the restricted space. 
Furthermore, the hopping elements depend on the elec- 
tronic configuration in c orbitals. In case of empty c 
orbitals, as e.g. in Sr2V02, the hopping elements in Eq. 
(|2.2p would be given by the bare tight binding element,— 
i.e., t a = t. On the contrary, the present situation with 
filled c orbitals in doped Yi^Ca^VOs perovskites re- 
sembles the case of doped manganites^ tA where the 
hopping elements between active e g orbitals are strongly 
renormalized by the tig spins S — 3/2 on both sites. In 
the manganites this follows from strong Hund's exchange 
coupling Jh between e g and t2 g electrons which stabilizes 
a HS S — 2 state at each Mn 3+ site.— Here one has in- 
stead a spin S = 1/2 of a c electron at each site which 
couples by Hund's exchange to the spin s = 1/2 of a sec- 
ond electron in the {a, 6} doublet, and a HS state with 
spin 5* = 1 arises. Therefore the constrainted hopping 
Hamiltonian T~it which follows from the symmetry of ti g 
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orbital states, 

#*= - * E («t« JCT +C^+H.c.) 

■ - 1 E + H - c -) - t E (°U*, + Hx -) ■ 

<ij>ll« 

(3.9) 

and contains the operators which act in the projected 
space, with electron number per site being either rii = 2 
or rii — 1. When two electrons are present (m — 2), they 
interact by the local exchange term 'HHund in Eq. (|3.2|) . 
which reads: 

'HHund = —2Jh ^(5ia + Sib) ' Si c . (3.10) 

i 

Altogether, Eqs. (pT9"J) and (|3~10l) define the double 
exchange model^ for strongly correlated {a, b} electrons 
interacting with localized spins s=l/2 of c electrons. In 
the effective model which follows from it and is described 
in Sec. [V] the effective hopping amplitude £y < t for the 
mobile {a, b} electron does depend on the directions of 
two electron spins in c orbitals along a bond (ij). We 
analyze the kinetic energy obtained in both AF phases 
below in Sec. [V] 

D. Orbital-charge interaction near Ca defects 

When an Y ion in Yi_a;Ca x V03 is replaced by a Ca im- 
purity, the lattice is disturbed and two interaction terms 
arise due to the presence of the impurity. The first of 
them is the Coulomb potential due to the Ca impurity, 
introduced before in Eq. (|2.10[) . while the second one is a 
crystal field term Hd considered here as the second last 
term in the spin-orbital t-J Hamiltonian Eq. (|3.2|) . The 
former term causes that a hole in the VO3 subsystem 
is confined to the immediate neighborhood of the Ca 2+ 
charge defect in the Y 3+ sublattice by the electrostatic 
potential Eq. (|2.10j) . as we have verified by numerical 
calculations reported in Sec. HU The latter term origi- 
nates from the quadrupolar component of electrostatic 
field generated by a charge defect at the V 3+ ions. Re- 
placing an Y 3+ ion by a Ca 2+ ion implies that an effective 
negative charge e is introduced at the Ca impurity, i.e., 
in the center of the cube shown in the inset of Fig. [6] 
It interacts with an electron in the {a, b} orbital doublet 
of the considered V 3+ ion, and the repulsive energy be- 
tween this ion and the Ca defect has to be minimized. 
This may be achieved by an optimal choice of the occu- 
pied tig orbital in the {a, b} subspace. 

Thus, the repulsive Coulomb potential of the Ca defect 
generates an orbital polarization at surrounding it V ions 
with a pronounced tendency toward electron occupation 
of one of the two linear combinations of the active {a, 6} 
orbitals at site i, 

|±)^i(al±6|)|0), (3.11) 




FIG. 6: (Color online) Energy levels of the t2 g orbital states 
7 = a,b,c for a representative V 3+ ion in YVO3 (left), and 
their modification in the vicinity of a Ca defect by the im- 
purity potential Him P (middle) and after adding the charge- 
orbital interaction Hd (right). Two orbitals a — 1 and 2 are 
occupied at V 3+ ions in the vicinity of the defect site D, as 
indicated by electron spins (arrows) which contribute to spin 
5 = 1. Occupied t-2 g orbitals a — 2 with intermediate energy 
are shown in the inset; there orbitals are linear combinations 
of {|o), |6)} orbitals |±)», defined in Eq. p. lip . Orbital phases 
and the c orbitals occupied at each site are not shown for clar- 
ity. In the doped case one of the occupied orbitals a — 2 is 
replaced by a hole in the vicinity of the defect D (see Fig. 1). 

being eigenstates of af Pauli matrix. A properly cho- 
sen linear combination, either or |— )j state at site 
i, maximizes the average distance between the interact- 
ing electronic charges at the V 3+ ion and at the Ca 2+ 
impurity, and minimizes the electrostatic interaction en- 
ergy. The orbital states which satisfy this condition and 
are favored at V 3+ ions are shown in the inset of Fig. [5] 
The lifting of orbital degeneracy will be described by the 
following orbital-defect polarization term acting in the 
{a, b} sector^ 

n D = d ]T vf , (3.12) 

where rf is defined in Eqs. (|3.7p . The summation over 
i € C\ in Eq. (|3.12p includes the V sites which belong to 
the cube C\ around the considered Ca site, and the sign 
factor Aj = ±1 (for D > 0) selects the proper orbital 
state which minimizes the charge-orbital interaction at 
each site i. This polarized orbital is labeled as a = 2 
in Fig. [6] It is occupied in V 3+ ions, but an extra hole 
introduced by doping does remove the electron from this 
topmost occupied orbital at a V 3+ ion. 

As an illustrative example of the expected conse- 
quences of the orbital-defect polarization interaction D, 
see Eq. (|3.12[) . we present the change of the orbital state 
for a single bond in Appendix [XJ One finds a first or- 
der transition from the orbital singlet with fluctuating 
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orbitals to a fully polarized orbital state at large D. We 
show below that although the considered bond interacts 
with its neighbors and orbital fluctuations are reduced 
from those in a single bond, the main features of this 
orbital transition survive also in a crystal. 



E. Superexchange for V 4+ -V 3+ bonds 

As already discussed in Sec. IIIQ doping by a hole 
transfers a V 3+ ion into a V 4+ ion and excitations to 
the UHB are then removed at the hole site, see Fig. 
[3j For the orbital background this implies a charge de- 
fect and one may consider the problem of hole propaga- 
tion, either in the orbital model^ 7 - or in the spin-orbital 
strong-coupling models Simultaneously, however, su- 
perexchange interactions change in a drastic way and 
may be treated in an analogous way as the superexchange 
between Mn 3+ and Mn 4+ ions in doped manganitesi 7 -^ 
Here the superexchange couples spin 5 = 1 of a V 3+ ion 
with spin s = 1/2 at the hole site, and may be derived by 
realizing that: (i) only charge excitations d\d 2 =± d?dj 
along the bond contribute and they depend on the orbital 
degree of freedom, (ii) only charge excitations to the LS 
(5 = 0) states contribute and can be treated in perturba- 
tion theory, and (iii) the excitation energies are given by 
Hund's exchange Jh as in the manganites^ There are 
two d 2 singlet states: (i) interorbital singlets with energy 
(U—Jh), excited in c|cj(a/6)j ;= c}(a/b)jcj processes,— 
and (ii) double occupancies of the c orbitals, excited in 
c}cj(a/b)j =± c 2 (a/b)}j processes. Due to the symme- 
try of the Coulomb interactions in the t2 g subspace, the 
latter excitations contribute via two eigenstates with en- 
ergies (U — Jh) and (U + 2Jh), respectively^ Therefore, 
taking the energy difference with respect to the energy 
(U — 3Jh) of the initial HS state at site j, one finds two 
excitation energies: 2Jh and 5Jh ■ 

Note that for the particular occupancy c 1 (a/6) 1 of ti g 
orbitals realized for V 3+ ions in YVO3, see Eq. (|2.3I) . 
the cf excited states are generated only in charge exci- 
tations along the bonds in the ab plane, while the for- 
mer interorbital c\(a/b)\ singlets may arise on the bonds 
along all three directions. Therefore, the superexchange 
between the hole site (V 4+ ion) and the neighboring un- 
doped (V 3+ ) sites is anisotropic. More details concerning 
the derivation of the superexchange between the s = 1/2 
spin of V 4+ ion and its V 3+ neighbors with 5 = 1 spins 
are presented in Appendix [5] 

While the superexchange in the ab plane depends in 
principle on the orbital {a, b} flavor at the d 2 site and on 
the bond direction, i.e., whether or not the hopping of 
the electron with either a or b orbital flavor is allowed (or 
not) along the considered bond (ij), see Appendix [Bj we 
present here the superexchange after averaging over the 
orbital configuration at the d 2 site. This simplification 
is well justified as we consider below (in Sec. [Vj the 
effective ID models along the c axis, where the spin and 
orbital configurations at the neighboring sites along the 



bonds parallel to either a or & axes are averaged out. 
Therefore, we introduce the exchange constants: 

t 2 ( t 2 2t 2 \ 

'-aw' Mew + ew)' (3 ' 13) 

and write the superexchange terms as follows, 
Hi = I c ^2 " 83 ~ 2) i 1 ~ n i) n 3 

+ lab J2 ( * ' 3j ~ \ ) _ • ( 3 - 14 ) 

(<3>||o6 ^ ' 

Here the spin operators s*j refer to s — 1/2 spin at the hole 
site), while Sj stands for 5=1 spins of the undoped V 3+ 
sites neighboring with the hole site, so the superexchange 
energy contributes for the bonds between V 3+ and V 4+ 
ions, i.e., when (1 —m)nj = 1. This result is used below 
(in Sec. |VJ) to derive and investigate ID orbital models 
which contain one doped hole. 

IV. ONE-DIMENSIONAL ORBITAL PHYSICS 
A. Motivation and calculation method 

In this section we analyze the perturbation of the spin- 
orbital structure due to the presence of defects, both for 
the G-AF and the C-AF phase. The main difference 
between the two phases arises from the magnetic order 
along the c axis, being either AF (FM) in the G-AF (G- 
AF) phase of Yi-^Ca^VOa. On the contrary, in the ab 
plane both phases share a common AF and AO struc- 
ture. Thus for our estimate of the difference of the free 
energies the orbital correlations along the c axis will be 
of particular importance. 

The microscopic model Hamiltonian given in Eq. (|3.2[) 
is too complex to treat simultaneously spin and orbital 
dynamics in a controlled approximation. Previous work 
has shown that larger 5=1 spins have weaker quantum 
effects and could be treated classically, in contrast to 
the orbital r = 1/2 pseudospins, which undergo strong 
orbital quantum fluctuations in the G-AF phase along the 
c axis, and play a crucial role to explain its anisotropic 
magnetic properties the temperature dependence of 
the optical spectral weights^ and the phase diagram of 
the -RVO3 perovskitesJ^ We emphasize that the bonds 
in ab planes are AF and thus similar in both magnetic 
phases. To capture the difference between the G-AF and 
G-AF phase we focus here on the ID orbital models along 
the c axis derived from the microscopic spin-orbital model 
Eq. (|3.2p . These orbital models are solved below for 
finite ID clusters coupled to the mean field (MF) terms 
which arise from the bonds in ab planes and capture the 
main difference between both magnetic phases, the G- 
AF and G-AF phase of Yi-^Ca^VOs. We analyze these 
models below in Sees. HVI andlVl Note that this approach 
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allows one to include the leading quantum fluctuations in 
the ground state; a similar cluster calculation was used 
recently to investigate the phase diagram of the Kugel- 
Khomskii model for a bilayeriS 

The orbital chains derived from the full spin-orbital 
model Eq. Q3.2j) are solved self-consistently using the MF 
terms arising from the interactions with the neighboring 
V 3+ ions along the bonds parallel to a and b axes. This 
simplification arises when bond orbital correlations are 
treated classically and replaced by their MF values. The 
orbital order parameter is then defined as 



X 



N 



(4.1) 



where X — G, G and Qc — (0, 7r,7r), Qq — (t, tt, tt) 
are vectors from the reciprocal space that correspond to 
the orbital alternation in the C-AO and G-AO phase, 
respectively. This approach is well justified here as the 
interactions in ab planes are Ising-like. Furthermore, it 
allows one to focus on the quantum fluctuations along 
the ID orbital chains parallel to the c axis, where both 
{a, b} orbital flavors are active, and on the role played by 
the orbital polarization term Eq. (|3.12l) . 



B. Orbital chain in the G-AF phase 

We consider first the undoped G-AF phase with clas- 
sical spin order, 



(Si ■ 



-1, 



(Si ■ Sj) a i) 



1 



(4.2) 



and C-AO order, stable at low temperature in undoped 
YVO3. Orbital interactions in the ab planes will be in- 
cluded in form of MF terms which contain contributions 
both from the superexchange oc J and from the Jahn- 
Teller-type orbital interactions oc V a b- These two inter- 
actions support each other, similar to the situation en- 
countered in LaMn03.— For the present G-AF phase with 
G-AO order the MF terms acting on every site of an or- 
bital chain along the c axis are the same and have alter- 
nating sign between two neighboring chains. We consider 
here a representative orbital chain with the G-AO order 
with the orbital order parameter, (t z )c > defined in 
Eq. (14.11) and stabilized by an effective field, 



h c = {2J V (r 1 +r 3 )+AV ab }(T z ) c . 



(4.3) 



The field originates from four bonds perpendicular to the 
chain and belonging to an ab plane and acts on the or- 
bital pseudospins {rf } within the chain, see below. The 
effective ID orbital model for the orbital chain along the 
c axis within the G-AF phase takes therefore the form, 



H, 



ID 



1 N 

- J (2ri + 1 - 3r 3 ) ^ n ■ r i+1 



N 



i=l 
2 



N 



(a) 
1< 



+> 



+> 



(b) 
1< 



(c) 



+> 



2 t 

1 

© Ia> 

4© l a > 
c © Ia> 



6 <)|a> 

G-AF 
C-AO 



D 



+> 



w 



2* 

3 ?' b> 

s 

C |a> 



O' 

... 

1 ^ 
1 

9 



w 



9 



6 

C-AF C-AF/G-AO 
G-AO +hole 



FIG. 7: (Color online) Schematic view of a section of an 
orbital chain parallel to the c axis. A Ca defect D in the 
vicinity of the chain is shown in the top part, close to the 
bond (12), and marked by dark filled circles; more distant 
nonequivalent sites are labeled as i = 3, • • • , 6. The orbital 
chains correspond to: (a) the G-AF phase with a electrons 
in I a) states (dark semicircles) corresponding to G-AO order 
along the chain (here sites i = 1 and i = 2 are equivalent); 
(b) the G-AF phase with G-AO order represented by alter- 
nating {|a), |6)} occupied orbital states along the chain (dark 
semicircles); (c) the G-AF phase with G-AO order of {a, b} 
orbitals as in (b) and a hole h at site i = 1 near the defect D. 
The orbitals on two top sites i = 1, 2 belong to the Vg cube 
around Ca defect D (see Fig. [6|, and are modified toward the 
|+) orbital states by the increasing orbital polarization inter- 
action D, see Eq. (|3. 12[) . The hole in (c) is confined to sites 
i = 1,2 by the trapping Coulomb potential Ed- Away from 
the defect D, orbital correlations in (a) are (almost) uniform 
and support AF spin coupling, while in (b) and (c) they gen- 
erate instead FM interactions that alternate between strong 
(s) and weak (w) exchange bonds. 



The sign selected in the last term oc he gives indeed 
a positive orbital order parameter (t z )c Eq. (|4.1|) 
when calculated self-consistently for the considered or- 
bital chain. The charge-orbital interactions oc D acts 
only at the bond (12) which belongs to a cube surround- 
ing a Ca site, i.e., at sites i = 1 and i = 2, and favors the 
orbitals shown in Fig. [6] A schematic view of the orbital 
chain in the G-AF phase is presented in Fig. [7ja). 

We emphasize that only finite orbital interaction V c > 
accompanied by the planar field oc he > stabilizes 
the observed G-AO order in the G-AF phase, while the 
superexchange term alone (at V c = and he — 0) would 
favor instead alternation of occupied {a, b} orbitals for 
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FIG. 8: (Color online) On-site orbital expectation values for 
the G-AF phase for increasing charge-orbital interaction D, 
see Eq. (|3.12[) : (a) orbital polarization (rf), and (b) or- 
bital order parameter (jf ). Labels of different curves indicate 
nonequivalent sites i = 1, 3, 4, 5 in a chain of L z — 8 sites, see 
Fig. Eta). Parameters: V ab = 0.2 J, V c = 0.7 J. 



finite Hund's exchange rj > 0, with (f^ • 7}+i) < 0. For 
the parameters of Fig. [5] one finds that orbital fluctua- 
tions along the orbital chain are quenched in the C-AO 
order and the FO order is classical with (t z )c = 0.5 in 
absence of charge-orbital interaction. This result is inde- 
pendent of the chain length L z . We show below how the 
classical C-AO order is modified by finite charge-orbital 
interaction D which we treat here as a free parameter. 

When a charge defect is introduced, the FO order in 
the chain is locally disturbed, and the occupied orbitals 
on the two equivalent sites i = 1 and i = 2 of the bond 
(12) gradually rotate with increasing D toward the linear 
combinations which minimize the charge-orbital interac- 
tion, see the inset of Fig. [6l One thus expects that a final 
state obtained for sufficiently large polarization interac- 
tion D has fully polarized orbitals, as shown schemati- 
cally in Fig. E{a). The FO order along the orbital chain 
parallel to the c axis remains undisturbed away from the 
Ca defect, while close to it the occupied orbitals change 
to |+) orbitals for the sign selected in Eq. (|4.4[) . 

The above scenario was confirmed by an exact diag- 
onalization of orbital chains of size L z — 8 sites which 
demonstrated local nature of the perturbation introduced 
by a Ca impurity. The orbital state near the charge de- 
fect changes gradually under increasing interaction D, as 
shown in Fig. [5] (We have verified that the length of 



L z — 8 sites is sufficient and longer chains of L z — 10 or 
12 sites give practically the same numerical results.) For 
the present case of classical C-AO order, the bond (12) 
decouples easily from the chain and the occupied orbitals 
at sites i = 1,2 are rotated, as shown by increasing ex- 
pectation value of (rf ), see Fig. [S^a). This local change 
near the defect is accompanied by weak negative polar- 
ization at the next nearest neighbors of the defect (site 
i = 3), and has practically no effect at more distant sites 
i = 4 and i = 5. The orbital moments near the defect, 
(rf) = (t|), are gradually suppressed only near the de- 
fect, while the (rf ) ones remain close to 0.5 away from it 
(for i > 2) [Fig. Hfb)]. Therefore, the AF spin interac- 
tions along the chain, supported by FO order, are only 
weakened between the first and second neighbor of the 
defect site, i.e., at the bond (23) in Fig. [TJa). Altogether, 
the numerical results demonstrate that the orbital state 
is disturbed only locally near a Ca defect and these local 
modifications of the orbital state do not destabilize the 
coexisting G-AF and C-AO order away from the defect. 



C. Orbital chain in the C-AF phase 

In order to derive the form of the orbital chain for the 
C-AF phase we use the following spin correlations in the 
classical state: 



(Si ■ Si+i) c — 1 , (Si ■ Sj) a b — —1 



(4.5) 



The FM spin correlation function along the c axis sup- 
presses then all low-spin contributions to the superex- 
change Eq. (|3 . 3[) and the orbitals participate in strong 
quantum fluctuations along the chain, induced by the 
SU(2)-symmetric interaction oc fi ■ Ti+i with a large pref- 
actor Jr\. The model Hamiltonian which describes the 
ID orbital chain in the C-AF phase shown in Fig. [T^b) 
takes a similar form to the one for the G-AF phase dis- 
cussed in Sec. IIVBI 
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The G-AO order is here stabilized again by the planar 
field due to the vanadium neighbors in the ab plane: 



h G = {2Ji 1 {r 1 +r- i )+W ab }(T z ) G . 



(4.7) 



which is here proportional to the orbital order parameter 
(t z )g, alternating along the chain, see Eq. (14.11) . How- 
ever, in contrast to the classical FO order considered in 
Sec. IIVBI the orbital order is here rather weak as it 
competes with the orbital fluctuations along the chainJ^ 
Thus the ground state is stabilized by strong {a, b} quan- 
tum orbital fluctuations, while the real orbital order pa- 
rameter (t z )g is reduced. Moreover, the orbital interac- 
tions cx V c along the chain favor C-AO order and are in 
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FIG. 9: (Color online) On-site orbital expectation values ob- 
tained for the C-AF phase as a function of increasing charge- 
orbital interaction D, see Eq. (I3.12|) . obtained for a chain of 
L z = 8 sites: (a) orbital polarization (rf), and (b) orbital 
order parameter (r.f). Dotted line in (a) shows orbital polar- 
ization (t£ ) for the free chain (at h a b = 0). Labels of different 
curves and parameters as in Fig. [8] 



conflict with the G-AO order in the present case. This 
reduces the orbital order parameter (t z )q further as the 
field he is considerably smaller than he, see Eq. (|4.3p . 

First, we performed a self-consistent calculation for a 
reference chain of L z = 8 sites in absence of Ca defects, 
i.e., with no polarization term (D = 0), which gives the 
orbital order parameter (t z )g — 0.296 for the present pa- 
rameters, see Fig. HJb). Next we considered the orbital 
chain given by Eq. (|4.6[) near the charge defect, with 
orbitals polarized near the defect (D > 0) and coupled 
to the neighboring sites in ab planes by the MF terms. 
The relevant G-AO order is shown in Fig. EJb). Due 
to the strong {a, b} orbital fluctuations along the chain, 
the orbital state changes here much slower with increas- 
ing charge-orbital interaction D than in the G-AF phase, 
see Fig. |Hl The increasing orbital moments (rf ) = (rf) 
near the charge defect gradually disturb the fluctuations 
along the chain and induce also finite (rf ) moments on 
more distant sites which alternate along the chain, see 
Fig. Oa). In contrast, the (r?) moments are first undis- 
turbed as long as the change in the orbital ground state 
is small, but next they are gradually suppressed at the 
nearest neighbor sites of the charge defect when D in- 
creases beyond D > J, cf. Appendix |XJ Away from the 
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FIG. 10: (Color online) Orbital correlation functions {■fi-fi+i) 
found for the orbital chain with L z = 8 along the c axis in the 
C-AF phase shown in Fig. [TJb): (a) self-consistent calculation 
with the MF term oc he Eq. (|4.3[) acting on the chain from 
the neighboring sites in the ab planes; (b) free orbital chain 
{he = 0). Parameters: V ab = 0.2J, V c = 0.7J, t = 6.25J. 



defect the weak orbital order is practically undisturbed 
due to the MF terms, while near the charge defect, i.e., 
on the bond (12), the orbital state is locally modified by 
the charge-orbital interaction for D > 2 J. 

Further evidence that the bond (12) close to the charge 
defect decouples from the orbital chain when the charge- 
orbital interaction is sufficiently large is given by the or- 
bital correlations (fi ■ Ti+i). For the present parameters 
one finds {fi ■ n+i) — —0.437 for the reference chain of 
L z = 8 sites in absence of charge defect {D — 0), see Fig. 
[TOT a). The orbital fluctuations are only weakly reduced 
by the G-AO order from the value of {fi ■ r^+i) — —0.452 
obtained for a free chain, see Fig. ITOTb). 

When charge-orbital interaction is introduced {D > 0), 
the orbital fluctuations near the defect, {f\ ■ fty, are first 
reduced by the perturbative positive term cx D 2 when 
D -C J, but soon (at D ~ J) their reduction becomes 
close to linear in D, and finally this bond decouples from 
the orbital chain (for D ~ 2 J) as the orbital fluctuations 
are locally suppressed by uniformly polarized orbitals, 
see Fig. [TOTa). Therefore, the correlation function {f\ ■ 
T2) changes sign at D ~ 2.15J and becomes positive for 
larger values of D, approaching the limit of fully polarized 
orbitals, with FO correlation near the defect, (ri • T2) = 
0.25. At the same time, the remaining correlations are 
only weakly modified and correspond approximately to 
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those obtained for a shorter orbital chain with open ends 
under the influence of the MF terms. The largest change 
is found for the (f^ • T3) correlation which approaches the 
classical value —0.25 when D > 4 J (not shown), but the 
remaining orbital correlations are only weakly influenced 
by the presence of the charge defect. 

The essential feature of the evolution of orbital cor- 
relations with increasing orbital-charge interaction D is 
the gradual decoupling of the bond (12), next to the 
charge defect, from the fluctuating orbital chain. This 
phenomenon is more pronounced for a free orbital chain, 
shown in Fig. ITOT b). In this case one finds a some- 
what surprising effect of induced dimerized correlations 
along the chain, as explained below. In addition to the 
perturbative regime of weak D < J, where the orbital 
correlations are modified by terms cx D 2 , one recognizes 
here two distinct regimes, separated by a critical value 
D c ~ 2.39J at which the ground state changes abruptly 
and the bond (12) polarized by the charge-orbital in- 
teraction decouples from the fluctuating chain, see Fig. 
[7jb). At D = D c a quantum transition takes place from 
a jointly fluctuating orbital chain to broken chain with 
its ends neighboring with a static orbital-polarized bond 
near the charge defect. We remark that this transition 
from fluctuating to polarized orbitals near the charge de- 
fect is similar to the one which takes place for a single 
bond, see Appendix [XJ Common features are that the 
transition is also discontinuous here and occurs to the 
state with the orbitals near the charge defect being al- 
most fully polarized, with (t\ ■ T2) > 0.17. 

The two regimes separated by the above transition are 
quite distinct. For D < D c the charge-orbital interac- 
tion disturbs orbital fluctuations along the chain as they 
become gradually suppressed on the bond (12), and trig- 
gers dimerized orbital correlations. Orbital fluctuations 
are then enhanced on the bonds neighboring with the 
disturbed bond, represented here by the bond (23), see 
Fig. [T^b). This perturbation generates alternation be- 
tween weaker and stronger orbital correlations along the 
chain. In contrast, for D > D c the orbitals freeze in the 
polarized state on the (12) bond, and thus the orbital 
fluctuations become restricted by the constraint imposed 
by this bond at both ends of the remaining chain. Thus 
the role of stronger and weaker fluctuating bonds along 
the broken chain is now reversed, and the fluctuations on 
the bond (23) become weak. 

One finds that the above quantum transition for a free 
chain changes to a crossover for an embedded chain, see 
Fig. HPT a). Also in this latter case it leads to dimer- 
ized orbital correlations along the orbital chain but this 
alternation is much weaker than for a free chain. This 
might suggest that charge-orbital interaction could be re- 
sponsible for dimerization along the orbital chains, and 
would imply dimerized FM interactions in the C-AF 
phase. When spin interactions are evaluated for the or- 
bital correlations of Fig. [TUJb), FM interactions are in- 
deed weak on bonds (23) and (45) (shown by dashed 
lines) and strong on bonds (34) and (56) (shown by long- 



dashed lines). However, the coupling with other orbital 
chains along the bonds in ab planes reduces considerably 
these dimerized interactions, see Fig. ITUT a). Neverthe- 
less, we introduced average orbital correlation functions 
for strong and weak bonds, (fj • Ti+i)° and (t| • 7\+i)iu 
as shown in Fig. [3(b), and found that they differ from 
each other. We shall use them below in Sec. IVI Bl and 
demonstrate that the dimerization is indeed predicted by 
the present t-J orbital model Eq. (|3.2[) . yet its mecha- 
nism is more subtle. Apart from the orbital polarization 
in the vicinity of a defect which we have discussed here, 
the orbital order is also disturbed by a hole introduced 
with each defect. 



V. DOPED HOLE IN ORBITAL CHAINS 
A. Orbital t-J model for the G-AF phase 

We first analyze the local changes of spin and orbital 
correlations introduced by a hole in the G-AF phase. 
There are two distinct processes how the hole acts on 
the spin-orbital structure, namely via: (i) superexchange 
between s=l/2 and S=l sites, and (ii) double exchange 
which involves hole motion in the LHB (with spins in the 
HS configuration). The double exchange hopping pro- 
cess, i.e., (|, 1) — > (1, |), is controlled by the spin orien- 
tation of the involved spins of electrons in xy orbitals. 
One should note that two DE processes, one forward 
and one backward, are indeed distinct from the superex- 
change process (|, 1) —> (0, |) — > which involves 
an S — intermediate state with an excitation energy 
oc J n. Interestingly this latter proportionality makes 
this kind of superexchange larger than the conventional 
superexchange cx J, i.e., between S = 1 spins, which re- 
sults from virtual excitations across the Mott-Hubbard 
gap cx U. The formal aspects of this kind of exchange 
were discussed in Sec. MI El 

Due to large binding energy between the hole and the 
Ca impurity, the hole is confined to one of the topmost 
occupied orbitals of the V 3+ cube shown in Fig. [HJ As 
discussed in SeclUthe motion of the the hole along the 
bond (12) parallel to the c axis is of particular relevance. 
A schematic picture is given in Fig. EKc), where one of 
the polarized orbitals in Fig. [7](a) is replaced by a hole. 
The AF order in the ab plane is stabilized mainly by 
the superexchange driven by the excitations in c orbitals 
which arc unaffected by the doped hole. This implies that 
the hole motion is confined to the (12) bond along the 
c axis, and the effective hopping element is determined 
here by the double exchange. 

As usual in the double exchange model^ the AF order 
is disturbed and spins cant in order to find a compromise 
between the loss of the magnetic energy with respect to 
the AF spin order, and the kinetic energy, which would 
be minimal for FM spins. In the present model with hole 
confinement near the charge defect, it suffices to analyze 
the double exchange mechanism on the bond (12) along 
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FIG. 11: (Color online) Schematic view of the hole hopping 
along the bond (12) belonging to the cube surrounding a Ca 
defect for the: (a) G-AF phase, and (b) C-AF phase. In 
the initial state \i) two spins s = 1/2 due to electrons which 
occupy |c) and \a) orbitals form a high-spin S = 1 state at d 2 
site i = l. The d 1 site i = 2 is occupied by the hole (filled 
circle) and has only a single s = 1/2 spin due to the electron 
occupying the |c) orbital. In the G-AF phase shown in (a) a 
hopping process of the electron with a orbital flavor, oc ta[a 2 , 
generates a finite canting angle 9 for both spins with respect 
to their original orientation (8 — 0) — which leads to a finite 
kinetic energy gain oc isin# by the transition to the final |/) 
high-spin S = 1 state. In contrast, full hopping amplitude t 
contributes along FM bonds in the G-AF phase, see (b). 



the c axis, see Fig. [TTJ supplemented by the magnetic 
energy contributions arising due to other bonds which 
start at either i = 1 or i = 2. In the initial state \i) 
shown in Fig. Illf a) the spin of an a = yz electron at site 
i = 1 is parallel to the one of the electron in c orbital. 
The hopping process to the final state |/) with the a 
electron moved to site i — 2 is possible after the spins 
at both sites are canted by angle 9 away from their AF 
order in the G-AF phase, see Fig. fTTTa). In this case the 
hopping amplitude is given by 



t(9) = tsmO. 



(5.1) 



LU 




0.4 0.6 
sin9 



FIG. 12: (Color online) Total energy per site £a calculated 
for the G-AF state with G-AO orbital order as function of tilt 
angle sin#, as obtained by exact diagonalization of an orbital 
chain of L z = 8 sites with D = 2J (solid line) and D — 
(dashed line). Parameters: V a t = 0.2 J, V c = 0.7 J, t = 6. 25 J. 



(relevant for the C-AO order): 

Ke(°) = -t{6)(a\a 2 + b\b 2 + R.c.) 

N N 



i=i 



N 



(5.2) 



Consequently, the orbital chain containing a hole in 
the G-AF phase is described by the orbital Hamiltonian 



As in Eq. (|3.2p . the creation operators {a^, tf} of spinlcss 
fermions act in the restricted space, and the hopping oc- 
curs along the (12) bond. The orbital chain depends on 
the orbital MF he, see Eq. (|4.3|) , and on the spin canting 
angle 9 via the hopping t(9), as given in Eq. (|5.ip . 

We performed exact diagonalization of the orbital 
chain Eq. (|5.2[) for representative parameters of a doped 
Yi^^Ca^VOs system. Thereby we investigated the to- 
tal energy including the magnetic energy £™ ag (6>) which 
follows from the bonds which are influenced by the spin 
canting at sites i — 1 and i = 2, see Appendix C. As ex- 
pected, one finds that the energy is lowered when the spin 
order at the hole site and in its neighborhood is locally 
disturbed and permits hopping with the reduced hopping 
element given by Eq. (|5.1I) . The kinetic energy oc t(9) is 
gained and part of the magnetic energy Eq. (|C1[) is lost 
when the spins cant and rotate away from their orienta- 
tion in the ideal G-AF phase. The optimal angle is found 
to be given by sin6> ~ 0.42 at D = 2 J and sin6> ~ 0.43 
at D = 0, see Fig. [T2j Thus, the dependence of 9q on 
the charge-orbital polarization term is surprisingly weak 
and we may consider the angle sin6*o — 0.42 obtained 
for D = 2 J as a representative value for the analysis of 
doped Yi-^Ca^VOs, see Sec. IVIBI We have found that 
already rather weak charge- orbital interaction D < J 
polarizes almost entirely the occupied orbital within the 
(12) bond which decouples from the remaining part of 
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FIG. 13: (Color online) Orbital order parameter (rf } in the C- 
AF phase around the charge defect neighboring with the bond 
(12) shown in Fig. [TJc), obtained for an embedded orbital 
chain with finite field ha for increasing orbital polarization 
interaction D, see Eq. (|3.12|l . Parameters: V a t = 0.2 J, V c = 
0.7J and t = 6.25J. 

the orbital chain. 

It is worth noting that the spin structure is robust and 
its local modification is moderate. This contradicts naive 
expectations based on rather small magnetic exchange 
constants between S = 1 spins, being only a fraction 
of J and determined by neutron experiments^ that the 
kinetic energy would dominate over magnetic one since 
t ^> J. Instead, we have found that the canting of the 
spin structure is only moderate for the realistic parame- 
ters, see Fig. [12] This surprising result follows from large 
AF superexchange interactions in the vicinity of the hole, 
I c ~ 1.92J and I a b — 4.1J, which are enhanced by small 
charge excitation energies oc Jh- Here we used Jh = 0.64 
eV for rj = Jh/U = 0.13. 13 Indeed, these exchange con- 
stants (I c and I a b) are larger by approximately one order 
of magnitude than the exchange constants J c and J a b 
between S — 1 spins. 



B. Orbital t-J model for the C-AF phase 

In contrast to the G-AF phase discussed before, the 
spin structure of the C-AF phase is already tuned to op- 
timize the double exchange, i.e., the hole hopping term 
is t near the charge defect, see Fig. [TlTb). An intriguing 
question here, however, is to what extent the spin struc- 
ture is affected indirectly via the perturbations of the 
orbital chain. After introducing a hole, an orbital chain 
with an even number of sites as considered in Sec. IIV CI 
becomes an open chain with an odd number of sites and 
an additional constraint — the orbital at one of its ends 
is polarized by the charge-orbital interaction oc D, and 
this orbital may interchange its position with the hole. 

The orbital Hamiltonian which describes a hole in the 
C-AF phase (with G-AO order) is easily derived from Eq. 



(SD, 

He = -t{a\a 2 + b\b 2 + H.c.) 

||c (z,i+l) ||c 

- D r- - h G J2(-^r: ■ (5.3) 

i=l,2 i 

It depends on the MF Hq defined in Eq. (|4.7[) which 
stabilizes the G-AO order and is relevant for the C-AF 
phase. Here we adopted the usual notation that the bond 
(12) belongs to the cube surrounding the charge defect, 
see Fig. [7f c) . As the lone electron at the undoped site is 
delocalized over the bond (12), the corresponding orbital 
is practically decoupled from the orbital chain along the 
c axis, and is easily polarized to minimize the charge- 
orbital interaction by a rather weak interaction, D <C 
Therefore, the order parameter (r^) vanishes at this bond, 

1. e., for i = 1 and i — 2, see Fig. [13J One finds as 
well that the orbital moments (rf) are undisturbed and 
alternate in the remaining part of the chain, as expected 
for the G-AO order. 

The most important consequence of the fragmentation 
of the orbital chain by the presence of a doped hole is 
the alternation of orbital correlations (t^ • Ti+i) shown 
in Fig. Q3] The hole due to the confinement by the 
defect potential moves predominantly on the bond (12) 
which leads to strong orbital fluctuations on sites 1 and 

2. These strong fluctuations suppress the correlations 
on the neighbor bond (23) (and on the equivalent bond 
(IL Z ) under periodic boundary conditions). 

Therefore, bonds with enhanced and reduced orbital 
fluctuations |(Tj • t;+i)| alternate along the chain and, 
in analogy to Sec. IIV CI we introduce the average or- 
bital correlation functions for strong and weak bonds: 
(t« • Ti + i)h s and {fi ■ Ti + i)hw- This alternation is most 
pronounced in a free chain, see Fig. 114( b). but is suffi- 
ciently robust to survive the coupling of the chain to its 
neighboring sites in the ab planes of the C-AF phase, see 
Fig. 1141 (a) . Therefore, we find here dimerization of the 
orbital interactions around a hole which is the driving 
force toward dimerized FM spin interactions in a doped 
C-AF phase, see Sec. IVIBI We also note that the alter- 
nation of orbital bond strength induced by the holes in 
the C-AF phase is reminiscent to the bond alternation 
induced by large charge-orbital interaction D, see Fig. 
Mb)- 



VI. SCENARIO FOR THE G-AF TO C-AF 
MAGNETIC TRANSITION 

A. Statistical averaging in doped phases 

To demonstrate that a phase transition from the G- 
AF to C-AF phase indeed occurs with increasing doping 
x in the Yi-^Ca^VOs system, we use here a statistical 
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FIG. 14: (Color online) Alternation of the bond orbital cor- 
relations {fi ■ n+i) as obtained for the C-AF phase in the 
vicinity of the charge defect, see Fig. [TJc) , as function of or- 
bital polarization interaction D Eq. (|3.12[) . The correlation 
functions which alternate between strong ((34) and (56), long- 
dashed lines) and weak ((23), solid lines, and (45), dashed 
lines) bonds, are obtained for an orbital chain with L z = 8 
sites: (a) embedded with finite he (V a b ~ 0.2 J), and (b) free 
with he = 0. Other parameters: V c = 0.7J and t — 6.25J. 



approach designed for the weakly doped regime. In the 
dilute limit we can neglect interactions among defects. 
Thus it suffices to analyze configurations of orbital chains 
of a fixed length with no more than one hole. A single 
calcium charge defect which occurs in Yi_ a; Ca a; V03 in- 
fluences the orbital order at eight vanadium ions which 
surround it in a cube shown in the inset of Fig. |BJ 

We start with the G-AF phase with C-AO order. This 
situation will be used as a reference state below and we 
shall investigate deviations from it by introducing charge 
defects. All orbitals in the orbital chains considered 
above are to some extent influenced by the presence of 
the charge defect, as we have seen in the previous Sees. 
IIVI and [Vj These modifications occur in a different way, 
i.e., depending on whether the chain under consideration 
is undoped (Sec. IIV[) or contains a single hole (Sec. |V|) . 
Among the direct neighbor chains of a defect, one chain 
contains a hole and the remaining three chains are un- 
doped. 

Next we want to determine the probability to find a 
chain segment of length L z that contains a hole. To 
this end we introduce the specific volume of a defect 



Vd = L 2 L Z . We assume here that the cubic lattice of 
V ions has lattice constant c = 1, such that the defect 
concentration is x — 1/L 2 L Z , and the total number of 
c-chains threading the box is I? . As we are interested in 
the dilute limit here, we choose for the box height L z = 8 
which is also the length of chains studied by numerical 
diagonalization. As L is in the range 2, 3, • • • ,L Z , the ac- 
cessible concentration range is x € (0.002,0.03). Smaller 
values for L z would allow us to study larger doping con- 
centration x. This may now be used to construct the 
statistical averaging over the bonds in a doped system. 
We introduce statistical weights or probabilities, 



w h = To 



1 



3 

T 2 



To > w 



L 2 -4 



(6.1) 



for finding a chain doped by a hole iwh) or an undoped 
chain next to a defect (w e ), respectively. Chains sep- 
arated from a cube with a defect site occur with the 
complementary probability wq = 1 — Wh — w e . These 
weights will be used below to: (i) estimate to what ex- 
tent dimerized interactions develop in the doped C-AF 
phase (in Sec. rvTBj) . (ii) obtain the energy of both AF 
phases which compete with each other, and (iii) inves- 
tigate whether a phase transition could occur from the 
G-AF to the G-AF phase in the low doping regime, see 
Sec. rVTCl 



B. 



Dimerization in orbital chains 



As we have shown in Sec. IV Bl a doped hole breaks the 
orbital chains, suppresses locally the orbital fluctuations 
and thus induces dimerized orbital correlations (see Fig. 
I15p . To some extent a similar dimerization occurs also in 
the undoped orbital chains that are direct neighbors of 
charge defects (see Fig. [TU)) . Although the dimerized or- 
bital correlations occur in the chain that contains a hole, 
we assume that the magnetic order and excitations will 
reflect an average alternation in the doped Yi^Ca^VOs 
crystal. Therefore, we average here the orbital corre- 
lations which occur on stronger and weaker fluctuating 
orbital bonds over the entire sample. 

The orbital correlation function for the undoped chain 
(tj • fi+i)^ = —0.437 is modified when doping x in- 
creases, and the fraction of bonds with stronger and 
weaker fluctuations gradually increases. Using the 
stronger and weaker orbital correlations calculated for 
the G-AF phase both for an undoped chain near the 
charge defect, (fj ■•7^+i)es(e«j) (Sec. IIV C[) . and for a chain 
containing a hole, (f, ■ ha(hm) (Sec. IV B[) . the aver- 
age orbital correlations for strong and weak bonds can 
be evaluated as follows: 

(Ti-T i+ l) s = (1 - W h - W e )(fi ■ T l+ l) (0) 

+ w e (fi ■ f l+1 ) es + w h (fi ■ f i+ i) hs , (6.2) 

(n-f i+ l) w = (1 - W h - W e )(fi ■ T l+ l) (0) 

+ w e (fi ■ f i+ i) ew + w h (n ■ f l+1 ) hw (6.3) 
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FIG. 15: (Color online) Dimerization in the C-AF phase in- 
duced by charge defects as function of the orbital polarization 
interaction D Eq. (|3.12|1 : (a) orbital correlations (fi-fi+i) on 
strong (i = 3, 5, solid lines) and weak bonds (i — 2, 4, dashed 
lines), and (b) magnetic exchange constants for strong (J cs , 
solid line) and weak (J cw , dashed line) FM bonds along the c 
axis. Dotted lines indicate the average values in both cases. 
Parameters: V ab = 0.2J, V c = 0.7 J, t = 6.25J, D = 3J. 



The difference between the stronger and weaker bonds 
increases with doping as shown in Fig. [Tola). At the 
same time the average orbital correlation function, 



Ti+l 



n+i) 



Ti+l) w ) , 



(6.4) 



increases somewhat with x, indicating overall reduction 
of orbital fluctuations with increasing doping. 

The FM exchange constants along the c axis also alter- 
nate when the orbital state in dimerized, see Fig. [T5fb) . 
They follow from the superexchange term in the spin- 
orbital model Eq. (|3.3[) by inserting average orbital corre- 
lations on strong (weak) bond into the formula for the ex- 
change constant 12 which follows from the superexchange 
along the c axis: 



J cw — 



r\r x 

- {n ■ T i+ i 



ri - rj( r i + r 3)] 
2 



:V r 3\ T i ■ T i+i 



)s}, (6.5) 



[ri - ri( r i + r 3)] 
2 



H+l/w 



:V r 3[Ti ■ t 1+1 



)»}.(6.6) 



Anisotropy between the exchange constants J cs and 



J cw , shown in Fig. fT5Tb). is caused by the alternating 
orbital correlations between strong and weak bonds, and 
increases with increasing doping x. The exchange con- 
stants may be thus expressed by the average exchange J a 
and the anisotropy 8 C > as follows: 



Jcs = Jo(l + S c ) , 

Jew Ja{\ $c ) • 



(6.7) 



The average FM interaction J a < is gradually weak- 
ened by reduced orbital fluctuations in a doped system. 
One finds J a = -0.245J and J a = -0.229J for x = 
and x = 0.02, respectively. In the present approach the 
anisotropy increases linearly with doping x and amounts 
to 5 C — 0.117 at x = 0.02. This anisotropy is much 
weaker than that found in the undoped YVO3 at T = 85 
K, where S c = 0.35. 1 — However, one should keep in mind 
that the above experimental result concerns the undoped 
compound (x — 0), where the present mechanism of 
dimerization is absent and only thermal spin fluctuations 
contribute. One expects that the present anisotropy of 
the FM exchange constants would be enhanced by the 
interchain coupling and by thermal fluctuations, both ef- 
fects not included in the present approach, and it would 
be then higher then in the undoped compound YVO3. 



C. Transition to dimerized C-AF phase 

The energy of the doped G-AF and C-AF phase was 
analyzed in a similar way — we evaluated the energy 
changes with respect to the undoped phases by weight- 
ing the terms arising for the orbital chains near the 
charge impurity, either with or without the hole, and 
we used thereby the weighting factors Wh and wq, see 
Eqs. (|6.ip . In the undoped YV0 3 (at x = 0) one finds 
the G-AF phase, but the energy of the C-AF is only 
larger by a rather small energy ~ 1 meV per site, as 
estimated before^ The magnetic excitations^ and op- 
tical experiments-! 2 suggest that J ~ 30 meV. We there- 
fore define the energy of the C-AF phase at x = as 
Ec(0) — 0.03 J and use it as a parameter below. 

The energies for both phases are given by: 



E G {x) 
E c {x) 



w e E% + w h E h G 



G 1 



w o E c {0) + w e El 



w h E h c . 



(6.8) 
(6.9) 



where E G 



£g{8o)- The energy increments E e x and E^, 
with X — G, C, follow from the analysis of the respec- 
tive orbital chains without a hole {E e x ) and with a single 
doped hole (E\). The magnetic contributions to the en- 
ergy and the energy difference originating from different 
orbital order in both undoped phases is already included 
in Ec(0). The magnetic order is only influenced locally 
near the charge defect in the chains containing a doped 
hole, and otherwise remains undisturbed, see Sec. IV Al 
Therefore we included also the magnetic term in Eq in 
Eq. dC5|. 
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x 

FIG. 16: (Color online) Energies of two competing magnetic 
phases as function of doping concentration x: G-AF phase 
(dashed line) and C-AF phase (solid line) obtained from 
Eqs. (|6.8[ ) and (|6.9|l . respectively. Parameters: V a b = 0.2 J, 
V c = 0.7J, t = 6.25J, E c {0) = 0.03J. 



By evaluating the energy contributions to Eg{x) and 
Ec{x) we derived the energies of both phases shown in 
Fig. (TBI One finds that the energy of the C-AF phase 
decreases faster than that of the G-AF phase when dop- 
ing increases. This results predominantly from the full 
contribution of the hopping t on the bond (12) occupied 
by the hole in the C-AF phase, while only a fraction of 
the kinetic energy t sin 6$ is released when it has to com- 
pete with robust AF order which hinders it in the G-AF 
phase. 



VII. DISCUSSION AND SUMMARY 

In this paper we addressed the microscopic reasons re- 
sponsible for the fast disappearance of the fragile G-AF 
phase with doping in Yi_ x Ca x V03 compounds. Our 
investigation of the electronic structure in the C-AF 
phase demonstrated that the undoped system is a multi- 
orbital Mott-Hubbard insulator and the upper Hubbard 
band consists of three subbands. The multiplet split- 
ting corresponds to excitations identified in the optical 
spectroscopy] 12 ' 15 i 57 Next we have considered doped sys- 
tems and we introduced a model for generic defects in 
the perovskite structure. The model is consistent with 
the experimentally observed variation of the optical spec- 
tra which show a distinct absorption deep in the Mott- 
Hubbard gap<^ Our interpretation of this maximum as 
originating from the vanadium states localized in the im- 
mediate neighborhood of defects turned out to be con- 
sistent with independent estimates of the binding energy 
of a hole close the Ca ion. This was further supported 
by the spectra obtained by analyzing the occupied and 
unoccupied vanadium states within the Hartree-Fock ap- 
proximation. In this way we derived large binding energy 
of the hole being ~ 1 eV, which guarantees that the hole 
is confined to a cube built by vanadium ions around the 
Ca defect. 



By reconsidering the local interaction parame- 
ters defining the multiband Hubbard model for the 
Yi-^Ca^VOs compounds and by analyzing the spectral 
functions obtained in the Hartree-Fock approach, we ar- 
rived at the conclusion that electron correlations are suf- 
ficiently strong to justify the description of doped ma- 
terials with the effective spin-orbital t-J model. This 
model was presented and analyzed for both competing 
AF phases: the G-AF phase stable at low doping, and 
the C-AF phase which takes over above the critical dop- 
ing x c ~ 0.02. The crucial observation is that the kinetic 
energy of the hole is controlled by the spin configura- 
tion, i.e., by a double exchange mechanism^ similar to 
that responsible for the ferromagnetic phase in doped 
manganites/2r— This makes the hopping processes in the 
two phases distinct along the c axis, where one phase has 
FM and the other one AF spin correlations. In the ab 
plane both phases have similar AF correlations, i.e., hole 
hopping is suppressed here by the spin order. Thus, the 
main difference arises from the hopping along the c axis 
which is not hindered in the C-AF phase and favors this 
phase. However, spins in the G-AF phase may adjust 
by tilting due to the double exchange and allow also for 
substantial gain of the kinetic energy. One might expect 
that the spins would cant to an almost FM alignment 
to favor the kinetic energy as t ^> J, where J/4 is the 
typical AF superexchange constant between 5 = 1 spins, 
which also determines the spin waves in the G-AF phase 
ofYV0 3 . 

Yet the double exchange mechanism, i.e., the canting 
of spins, is not controlled by J but by the exchange in- 
teractions I c and I a b around the hole site. These inter- 
actions couple the s — 1/2 spin at a V 4+ ion with the 
neighboring S = 1 spins — they are much larger than 
the energy scale J as these processes arise from exchange 
interactions violating Hund's exchange Jjj but not cre- 
ating d 3 configurations which cost the energy U, see Eqs. 
(|3 . 13[) . Thus, the canting of spins is opposed by rather 
strong AF interactions and the kinetic energy gain in the 
G-AF phase turns out to be only a fraction of the one 
obtained from the FM bonds in the C-AF phased As 
a result, the G-AF phase is stable in a window of low 
Ca doping, approximately for < x < 0.02. In this way 
we have identified the leading mechanism stabilizing the 
C-AF phase in doped Yi-^Ca^VOs which stems from 
spin-orbital physics: free hole hopping oc t on the bonds 
in the vicinity of the Ca defects accompanied by orbital 
fluctuations along the remaining FM bonds along the c 
axis. 

We have analyzed the differences in spin and orbital 
correlations around the defect states between the two AF 
phases: the G-AF and the C-AF phase. In the G-AF 
phase, stable at low doping, no orbital fluctuations can 
occur along the c axis due to the static nature of the 
C-AO order. It is for this reason that the orbital cor- 
relations are rather easily modified in the neighborhood 
of the Ca defect and the occupied orbital states follow 
the orbital polarization interaction imposed by the de- 
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feet. On the contrary, in the C-AF phase strong orbital 
fluctuations oppose the orbital polarization. Thus the 
orbitals are prevented to rotate toward the orientations 
preferred by the charge-orbital polarization interactions, 
as shown in Fig. [TOl 

As an important byproduct of the present study we 
have discovered that dimcrization of orbital correlations 
along the c axis is a natural consequence of doping and 
may be considered as induced by hole defects in the or- 
bital chains. In fact, a doped hole breaks the orbital chain 
and generates alternating orbital bond correlations. Due 
to the spin-orbital coupling, this alternation in the or- 
bital subsystem induces the alternation of stronger and 
weaker effective FM exchange interactions along the c 
axis. Therefore, the performed analysis leads to a pre- 
diction that dimerization of the FM exchange interac- 
tions should be enhanced by doping. We emphasize that 
this mechanism of dimerization is distinct from ther- 
mal fluctuations that are necessary to stabilize dimerized 
spin and orbital interactions in the undoped spin-orbital 
chain^ We expect that both mechanisms: (i) orbital- 
Peierls dimerization^ and (ii) defect induced dimeriza- 
tion of orbital correlations analyzed here, will support 
each other in doped Yi_ x Ca x V03. 

Summarizing, we have shown that a phase transition to 
the C-AF phase can be explained by the double exchange 
mechanism using the spin-orbital physics of the doped 
-RVO3 vanadates, as the kinetic energy of doped holes is 
much lower in the C-AF phase and this energy gain com- 
pensates the loss of the magnetic superexchange energy 
already at rather low doping x ~ 0.02. We have shown 
that the experimentally observed magnetic transition 33 
to the C-AF phase may be reproduced within the pre- 
sented microscopic model using the parameters consis- 
tent with other experimental data. 

Another challenging problem, not addressed here, is 
a surprising stability of the C-AF phase under dop- 
ing in Lai-^Sr^VOs compounds. The present calcu- 
lations have shown that large binding energy confines 
doped holes to the immediate neighborhood of Ca ions in 
Yi-zCa^VOa. If a similar hole confinement takes place 
also in Lai-^Sr^VOs, it could explain the insulating state 
found in these compounds in a broad range of doping 
< x < 0.18. 31 Although the theoretical explanation 
of the evolution of electronic and magnetic properties of 
Lai-^Sr^VOs has still to be constructed, we note that 
large binding energy between holes and defect states is 
consistent with these observations. Therefore, we suggest 
that it is a necessary ingredient of any realistic theoretical 
approach to the doped vanadium perovskites. 
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Appendix A: Orbital polarization transition 

In order to get a qualitative insight into the mechanism 
of local suppression of orbital fluctuations near the im- 
purity in the C-AF phase we consider a quantum transi- 
tion to the orbital polarization on a single bond (between 
sites i = l and i = 2 in Fig. [7]) from the singlet state 
fi-Tj — —3/4 to the fully polarized state with (rf ) = 1/2. 



The orbital Hamiltonian obtained from Eq. (|3.2p for the 
present toy orbital model with a bond (12) along the c 



axis is 



H = Jr% [tx ■ f 2 



D{rf 



(Al) 



We have chosen here the polarization interaction — D 
which favors (rf ) =1/2, but an equivalent result is ob- 
tained for +D. The Hamiltonian Eq. (|A1|) is easily di- 
agonalized using the basis {|w)}: 



|1) 
|2) 

|3> 
|4> 



44lo) 
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b\bl\0) 
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|0), 



|0), 



(A2) 



where |0) stands for the vacuum. One finds that the 
Hamiltonian matrix H nm = (n\H\m) takes the following 
form, 



V 








-7^ 

— j= D 










(A3) 



SJrJ 



and the triplet components {|1), |2), |3)} are coupled by 
the orbital polarization term cx D. The eigenvalues are: 



Ai 

^2,3 

A 4 



-Jr x ±D, 
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(A4) 
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As expected, the splitting between the singlet |4) and the 
triplet {|1), |2), |3)} state is Jr\ at D — 0. It decreases 
linearly with increasing D and one finds a quantum tran- 
sition at the critical value of polarization interaction, 



Jn. 



(A7) 
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This transition is first order and occurs as a level cross- 
ing between the singlet and fully polarized triplet com- 
ponent with energy A3 = \Jt\ — D. Above the transition 
(for D > Jr\) the orbital state is fully polarized by the 
charge-orbital interaction and (rf) = 1/2 for i = 1,2. 
In this orbital state fluctuations present in the orbital 
singlet are suppressed and triplet correlations take over, 
(n • f 2 ) = 1/4. 

The described orbital transition modifies also the mag- 
netic state in the coupled spin-orbital system. While the 
fluctuating orbitals in the singlet state support FM spin 
correlations, such a polarized orbital state at D > D c 
supports instead AF spin correlations along the bond in 
the spin-orbital model Eq. (|3.3|) . 




Appendix B: Superexchange for the d 1 — d 2 bond 

Here we present the derivation of the superexchange 
between the V 4+ ion generated by a doped hole and 
its V 3+ neighbors. When a hole is doped at a vana- 
dium ion in YVO3, the resulting spin electronic configu- 
ration is c\ (xy\), corresponding to S = 1/2 spin. The 
superexchange interaction follows from an interchange 
of charge between two sites in the excitation process, 
(c 1 ) i (c 1 (a/b) 1 )j ^ (c 1 (a/b) 1 ) i (c 1 )j, and involves only 
Hund's exchange Jjj due to the intermediate low-spin 
(S = 0) excited state (c 1 (a/b) 1 )i. Note that these ex- 
citations contribute to the AF superexchange, while the 
charge transitions between two ions in FM configuration 
do not involve any excitation energy and are treated ex- 
plicitly by the hopping term Eq. (|3.9p . similar as in doped 
manganites.— 

The actual occupancy c\ and cj(a/6)j of V 4+ and V 3+ 
ions is responsible for different contributions to the AF 
superexchange between the bonds along the c axis and 
in the ab planes. Consider first a bond (ij) \\ c. The 
excitations occur here solely by the hopping of an a (or 
b) electron to the neighboring site occupied by the hole 
and back, as shown in Figs. UTT a) andfTWb). The excited 
state, either c\a\ or cjb}, has to be next projected on the 
low-spin (S = 0) state, with the excitation energy of 2 Jh ■ 
As usual, the final state has the same charge distribution 
as the initial one, and the spin configuration is either the 
same as in Fig. H7T c)]. or the z-th components of spins 
at both sites have been changed by one, see Fig. UTT d). 
One finds 



4J ff 



3 2 



x {n»(l - rij) + rij(l — rii)} , 



(Bl) 



where rij is the number of electrons in the {a, b} orbital 
doublet, see Eq. (pT4"|) . 

The charge transitions which contribute to the su- 
perexchange along the bonds in ab plane, (ij) £ ab, have 
a richer structure as both ti g electrons at a given V 3+ 
ion may be allowed to hop to the hole site. While the c 



FIG. 17: (Color online) Artist's view of the virtual charge ex- 
citations d\d 2 — > d 2 dj —> djd 2 along a bond (ij) € C parallel 
to the c axis which contribute to the AF superexchange be- 
tween V 4+ and V 3+ ions in cubic vanadates. After a hopping 
of a electron from state (a) to the excited state (b), the exci- 
tation energy £1 = 2Jh arises. The excitation (b) may decay 
in two ways: (c) either the initial configuration is restored, or 
(d) spin flips take place at both sites. 



electron hopping is allowed along each bond, the second 
electron is either allowed to hop or not, depending on 
its flavor; we show in Fig. [T5]a case with both electrons 
contributing to the superexchange. First, the electron in 
the degenerate {a, b} orbitals may hop to the hole site, 
see Fig. [T8lb). This process is similar to the one for the 
bond along the c axis, but is allowed only for half of the 
bonds, depending on whether the hopping for the occu- 
pied orbital flavor is allowed or not (here we consider an 
a electron with the hopping allowed along the b axis) . It 
leads to two final states shown in Figs. \TE[c) andfTSTd). 
As a new feature one finds in addition the transitions by 
both c electrons which create a double occupancy in c 
orbital on the undoped site, cj, and leave behind the a 
(or b) electron [Fig. IToTe)]. This state has be projected 
onto the t^ g eigenstates with energies 2Jh and 5Jh , see 
Fig. 1 of Ref. 0. The final states, shown in Figs. USTf) 
and [PSf g). have again the same z-th spin states as the 
initial state, or the spins are flipped. One finds thus the 
superexchange, 



(«>||a V 

x {nib(l — rij) + rijb(l - n,)} 

+ BsE(*-4-f 

x {n ia (l - rij) + rij a (l - rii)} 

(i])\\ab 

x {rii(l - rij) + rij(l — m)} 



(B2) 
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=> 



=> 



(a) 



U Vt 

(b) => 



=> 



(e) => 



Ay Vt 

c) 



ty VI 

(d) 



iy Vt 

(f) 



V V± 

(g) 



d l —d 2 bond in the G-AF phase is obtained by minimizing 
the orbital problem Eq. (|5.2p together with the magnetic 
energy £™ ag (6l) which follows from the bonds which are 
influenced by the spin canting at sites i — 1 and i = 2. 
This energy consists of several contributions due to the 
bonds along the c axis and in the ab planes. When the 
considered site, either i = 1 or i = 2, is undoped, this 
energy is determined by the superexchange Eq. (|3.3[) 
between two S = 1 spins, while for the hole site it follows 
from the superexchange Eq. p,14|) between the spin s — 
1/2 at the hole site and its S = 1 neighbor spin. Using 
the MF approximation for the superexchange terms one 
finds: 



FIG. 18: (Color online) Artist's view of the virtual charge 
excitations d\d 2 — } dfdj — > d\d 2 along a bond (ij) € C along 
the b axis which contribute to the AF superexchange (a) be- 
tween V 4+ and V 3+ ions in cubic vanadates. When a elec- 
tron exchanges with a hole, the excited state (b) with energy 
£i = 2Jh arises — it leads to two ground state configura- 
tions, either (c) without or (d) with spin flip. The hopping 
of c electron creates a double occupancy (e) with the exci- 
tation energy ei = 5Jh — it gives again two ground state 
configurations, either (f) without or (g) with spin flip. 



The first two terms contribute only when the electron in 
the {a, b} doublet is allowed to hop along the bond (ij) 
in the ab plane, while the last term arises from the cf 
double occupancies and has no orbital dependence. 

Note that after the charge excitation the same orbital 
configuration has to be reached in the final state. There- 
fore, neither in this case (Fig. [T8|) . nor for the bond 
along the c axis considered in Fig. [T71 orbital fluctua- 
tions are present. The orbital dependence occurs in the 
interactions derived for the ab plane Eq. (IB2|) . but for 
the calculations for ID orbital chains along the c axis, 
performed in this paper, it suffices to average over the 
orbital configuration on the sites around the hole. Using 
the constraint Eq. (|2.3p one finds: 



(<:l)\\>-- 



rij - n%) 



(ij)\\ab x 



2t 2 
5Jh 



E (s i -S j - 1 -)(2 



rij - m) (B3) 



(ij)\\ab 

This result is used in Sec. [V] to investigate ID orbital 
chains which contain one doped hole. 



i=l,2 
i=l,2 (ij) || ab 

+ J?{(1 - n 2 )(S z 1 S z N ) + (1 - ni ){SIS§)} 
+ W-nz) E 

+ W-ni) E (W*>. (CI) 

(2j>||o6 



The hole disturbs the G-AF order locally, so the other 
bonds are only weakly influenced and one may evaluate 
the correlation functions in Eq. (|Cip using the classical 
spin order in this phase, see Eqs. (|4.2p . When the hole is 
at site i = 1, i.e., in the |/) state of Fig. [TTT a) (the other 
|i) configuration with a hole at site i = 2 is equivalent), 
we have used: 



(slSI) + (sfS 



A'/ 



i cos(26>) - i cos 6», (C2) 



(SISI) 



1 



cos 9 . 



2 

■ cos 6 . 



(C3) 
(C4) 



Finally, the ground state of the orbital chain containing 
one hole in the G-AF phase may be found by minimizing 
the energy obtained from the ID orbital chain Eq. f|5.2J> . 
including the correction of the magnetic energy £^ as (9) 



Eq. 4CT 



Appendix C: Energy evaluation in the double 
exchange model 



£ G (9) = {U%{6)) 



£^{9). 



(C5) 



Here we present the technical details of the calculations 
performed in Sec. IV Al The optimal canting angle for the 
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